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Abstract 

The geometric framework for the Hamilton-Jacobi theory developed in [HHIIES] is ex¬ 
tended for multisymplectic first-order classical field theories. The Hamilton-Jacobi problem 
is stated for the Lagrangian and the Hamiltonian formalisms of these theories as a particular 
case of a more general problem, and the classical Hamilton-Jacobi equation for field theories 
is recovered from this geometrical setting. Particular and complete solutions to these prob¬ 
lems are defined and characterized in several equivalent ways in both formalisms, and the 
equivalence between them is proved. The use of distributions in jet bundles that represent 
the solutions to the field equations is the fundamental tool in this formulation. Some ex¬ 
amples are analyzed and, in particular, the Hamilton-Jacobi equation for non-autonomous 
mechanical systems is obtained as a special case of our results. 
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1 Introduction 

The Hamilton-Jacobi theory, as it is classically presented in the textbooks and works on analyti¬ 
cal mechanics, is a way to integrate Hamilton equations (that is, a system of first-order ordinary 
differential equations), which consists in giving an appropriate canonical transformation leading 
the system to equilibrium [31331155]. This transformation is constructed from its generating 
function which, in this method, is obtained as the solution to a partial differential equation: 
the so-called Hamilton-Jacobi equation. This method is based on a famous contribution from 
Hamilton on geometric optics, where he showed that the propagation of wavefronts is charac¬ 
terized by a function (the characteristic function) which is the solution to a first-order partial 
differential equation called eikonal equation, which is related to the Hamilton-Jacobi equation. 
Thus, from a purely mathematical perspective, the Hamilton-Jacobi theory is a very important 
example showing the deep connection between systems of first-order ordinary differential equa¬ 
tions and first-order partial differential equations [54] • The Hamilton-Jacobi equation appears 
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also when short-wave approximations for the solutions of wave-type (hyperbolic) equations are 
considered. In this way, from a physical point of view, being a classical equation, it is also 
very close to the Schrodinger equation of quantum mechanics, since from a complete solution 
to the Hamilton-Jacobi equation, we are able to reconstruct an approximate solution to the 
Schrodinger equation [27( I48j and thus it allows to establish an approach within classical theory 
of the notions of wave function and state in quantum theory. 

For all these reasons, Hamilton-Jacobi theory is a matter of continuous interest, and it was 
studied and generalized also in other classical ambients; in particular, for constrained systems 
arising from singular Lagrangians (gauge theories) [2T] or also for higher-order dynamics |20j . 

Furthermore, in the last decades, great efforts have been done in understanding physical 
systems from a geometric perspective. Concerning to geometric mechanics, the intrinsic formu¬ 
lation of Hamilton-Jacobi equation is also clear and can be found in [U [46l I47j . In addition, in 
[14| a generic geometric framework for the Hamilton-Jacobi theory was formulated both in the 
Lagrangian and the Hamiltonian formalisms, for autonomous and non-autonomous mechanics, 
recovering the usual Hamilton-Jacobi equation as a special case in this generalized framework. 
In particular, it is shown that the existence of constants of motion helps to solve the Hamilton- 
Jacobi problem, which can be regarded as a way to describe the dynamics on the phase space 
of the system in terms of a family of vector fields on a submanifold of it. The basic ideas of this 
generalization of the Hamilton-Jacobi formalism are similar to those outlined in [35]. 

These geometric frameworks have been used by other authors to develop the Hamilton-Jacobi 
theory in many different situations in mechanics. For instance, the statement and applications 
of the Hamilton-Jacobi method for non-holonomic and holonomic mechanical systems is done 
in HKIIlllollMlIlIlllol, the geometric treatment of the theory for dynamical systems described 
by singular Lagrangians is analyzed in [371 H2l H3] , the application to control theory is given in 
PEUiig, and the generalization for higher-order dynamical systems is established in |18[ [19] . 
Moreover, the Hamilton-Jacobi theory has been extended for mechanical systems which are 
described using more general geometrical frameworks, such as Lie algebroids m i36], almost- 
Poisson manifolds [H], and fiber bundles in general and the relationship between the 
Hamilton-Jacobi equation and some other geometric structures in mechanics are analyzed in 
mm- Finally, the geometric discretization of the Hamilton-Jacobi equation is also considered 
in [5115I]. 

The extension of the Hamilton-Jacobi theory for first-order classical field theories has been 
developed for different covariant formulations (/c-symplectic and fc-cosymplectic) in the Hamil¬ 
tonian formalism mms] and also for the non-covariant Hamiltonian formulation (Cauchy data 
space) [To]. A first quick approach to state the Hamilton-Jacobi equation for the most general 
framework (the multisymplectic one) was outlined in [39], also in the Hamiltonian formalism. 
Furthermore, using a different approach involving connections, the theory has been generalized 
to higher-order field theories [59] and also for partial differential equations in general [60] . 

The aim of this paper is to complete these previous developments; that is, to use the guide¬ 
lines stated in the aforementioned references on the Hamilton-Jacobi theory in geometric me¬ 
chanics in order to give a complete description of this theory for the multisymplectic formalism 
of first-order classical field theories, both in the Lagrangian and the Hamiltonian formalisms, 
and showing the equivalence between them. Our standpoint is m and, in particular, some 
of our results are a development of the ideas pointed out in this reference. As a fundamental 
difference with these previous works, we consider the sections which are solutions of the field 
equations as integral sections of integrable distributions in the corresponding phase spaces (jet 
bundles and bundles of forms) where the equations are defined, and we represent these distribu¬ 
tions by means of (classes of) multivector fields in general [251126] . This allows us to adapt the 
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geometric models for the Hamilton-Jacobi problem in mechanics given in DU HU to the present 
case. 

The paper is organized as follows: Section [2] is a short review on multisymplectic geometry, 
jet bundles, and multivector fields and their relation with integrable distributions, which is given 
in order to do the paper self-contained. The main results are presented in Sections [3] and 01 
where first the generalized Hamilton-Jacobi problem, and later the standard Hamilton-Jacobi 
problem are stated in the Lagrangian and Hamiltonian formalisms, and the equivalence between 
both formalisms is analyzed. In these Sections, the particular and complete solutions of the 
Hamilton-Jacobi equations are introduced and interpreted geometrically. In Section [5l some 
examples are studied; in particular, non-autonomous dynamical systems as the particular case 
of a field theory with 1-dimensional base manifold, quadratic Lagrangian densities, and the 
problem of minimal surfaces in dimension three. Finally, the conclusions and further research 
are presented in Section [6l where the comparison and differences between our model and the 
aforementioned previous works are also discussed. 

All the manifolds are real, second countable and C°°. The maps and the structures are 
assumed to be C°°. Sum over crossed repeated indices is understood. 


2 Geometrical background 

2.1 Multisymplectic geometry 

In this section we give a short review on multisymplectic geometry and some particular sub¬ 
manifolds of a multisymplectic manifold (see [II1II21[22] for details). 

Let M be an m-dimensional smooth manifold. A multisymplectic k-form in M is a closed 
fe-form a; € Q^{M) which, in addition, is 1-nondegenerate, that is, for every p G M, i{Xp)ujp = 0 
if, and only if, Xp = 0, where Xp € TpM. If oj is closed and 1-degenerate, it is called a 
premultisymplectic k-form. A manifold endowed with a (pre)multisymplectic form is called a 
(pre)multisymplectic manifold of order k. 

Observe that a necessary condition for a fe-form to be 1-nondegenerate is 1 < A: ^ dim M. 

Given a symplectic manifold, we have a natural definition of “orthogonality” in terms of the 
symplectic form. This definition can be generalized to multisymplectic manifolds, bearing in 
mind that there are several levels of orthogonality to be considered. 

Definition 1. Let {M,oj) he a multisymplectic manifold of order k, and F C TM a vector 
subbundle. The Th orthogonal complement of F, with 1 ^ I < k is the subbundle C TM 
defined as 

= {{p,Up) G TM I u}p{up,vi,... ,vi) = 0 for every {p,Vi) G F} . 

Definition 2. A subbundle F C TM is called Lisotropic if F C F-^’\ Z-coisotropic if F^'^ C F, 
and Z-Lagrangian if F = F-^’\ for 1 ^ Z < Zc. 

Bearing in mind this last Definition, one can generalize the concepts of Z-isotropic, Z-coisotropic 
and Z-Lagrangian subbundles to immersed submanifolds as follows. 

Definition 3. Let {M,u}) be a multisymplectic manifold of order k, and N ^ M a submanifold 
with canonical embedding i: N ^ M. Let us consider the subbundle Ti{TN) C TM. Then, 
N is a Z-isotropic (immersed) submanifold (resp., Z-coisotropic submanifold, Z-Lagrangian sub- 
manifoldj ifTi{TN) is a I-isotropic (resp., I-coisotropic, l-Lagrangian) subbundle. 
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Finally, one has the following characterization of isotropic submanifolds of maximum order. 
Lemma 1. A submanifold i: N ^ M is {k — l)-isotropic if, and only if, i*oj = 0. 

2.2 First-order jet bundles 

In this section we give a short review on jet bundles: definition, some canonical structures and 
the concept of “dual bundle” (see 1571 for details). 


Definition and local coordinates. Prolongation of sections. Holonomy. 


Let M be an orientable m-dimensional smooth manifold with fixed volume form rj E f?™'(M), 
and let E M be a bundle with dimF = m + n. The first-order jet bundle of the projection 
TT, J^TT, is the manifold of the 1-jets of local sections (j) E r(7r); that is, equivalence classes of 
local sections of vr by the relation of equality on every first-order partial derivative. A point in 
J^vr is denoted by jl,4>, where x E M and f E r(7r) is a representative of the equivalence class. 
The manifold J^tt is endowed with the following natural projections 

TT^: J^vr —E vf^ = vr o TT^: J^vr —> M 

flf I—^ (j){x) ’ jlcj) I—^ X 


The hbers (tt^) ^(m) C J^tt, with u € E, are denoted Jjvr. 


Local coordinates in J^vr are introduced as follows: let (x*), 1 ^ i ^ m, be local coordinates 
in M such that r] = d^x = dx^ A ... A dx”*, and (x®,u“), 1 ^ a ^ n, local coordinates in 
E adapted to the bundle structure. Let (p E r(7r) be a section with coordinate expression 
0(x*) = (x*, (/>“(x*)). Then, local coordinates in J^tt are (x®, tt", uf), with 1 ^ ^ m and 

1 ^ a ^ n, where 


u 


d(f)^ 

9x® 


Using these coordinates, the local expressions of the natural projections are 


7r^(x®,u",uf) = (x®,u") ; 7r^(x®, nf) = (x®). 


\i (j) & r(7r) is a section, we denote the prolongation of cp to J^vr by j^cp E r(7r^). In the 
natural coordinates of J^n, if (p{x‘^) = (x®, i;f)“(x®)), the prolongation of cp is given by 

Definition 4. A section ip E r(7f^) is holonomic o ^p) = ip, that is, if there exists a 

section cp = o £ r(7r) such that ip is the prolongation of cp to J^tt. 

fn natural coordinates, if £ r(7r^) is given by 'ipiyJ) = {x^,'ip°‘, ipf), then the condition for 
Ip to be holonomic gives the system of partial differential equations 

dtp'^ 

C = l^i^rn, l^a^n, (1) 

An alternative characterization of holonomic sections is given in terms of the canonical 
structure form of J^vr. 
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Definition 5. The canonical structure form of J^n is the 1-form 6 in J^tt with values in V{tt) 
defined by 

where v € and is the vertical differential of (j) at (j){x) G E, and is defined as the 

^<p{x)E such that = T^(^)Id - T<^(a;)(<^ o tt). 

Proposition 1. A section G r(7f^) is holonomic if, and only if, 'fi*9 = 0. 


Dual bundles 

Let us consider the dual space of J^vr as an affine bundle over E, which is the set of affine maps 
from jfiT to (A™'T*M)jr(n)) with u € E, that is, the set 

U Aff(J>,(A-(T*M)),(,)). 

u&E 

From m we know that this set is a manifold diffeomorphic to the smooth vector bundle of 
TT-semibasic m-forms over E, Af'{T*E). This bundle is called the extended dual jet bundle of 
TT, and we have the following canonical projections 

7rE:Af-(T*E) —^ E _ 7rE:A^(T*F;) —^ M 

{u,uju) I —u ’ {u,Uu) I— x{u) 

Since Alf(T*E) is a bundle of forms, it is endowed with some canonical forms. First, the 
Liouville m-form, or tautological m-form, is the form 0 G n^{A^(T*E)) defined by 

Q{aj){Xi,.. .,Xm) = a;(T^7rs(Xi),... ,T^7rE(Xm)), 

where uj G Alf{T*E), and Xi, ..., X^ G Ttj(A™(T*Fl)). As usual, this form satisfies the property 
^*0 = ^ for every ^ G From this, the Liouville (m-\-l)-form, or canonical multisymplectic 

{m + l)-form, is the form Ll = —d0 G 17™'+^(A™(T*£')). 

Local coordinates in Af{T*E) are constructed as follows: let (x*) be local coordinates in 
M, and coordinates in E adapted to the bundle structure. Then, local coordinates in 

X!f{T*E) are (x*, where 1 ^ i ^ m, 1 ^ a ^ n. In these coordinates, the canonical 

projections have the following local expressions 

TTEixfu°‘,p,pl^) = {xfu^) ; 7rE(x*,u",p,p^) = (x*). 

On the other hand, the Liouville m and (m + l)-forms have the following local expressions 

0 = pd”^x + pidu^ A d^-^Xi ; Ll =-dp A d™x - dpi A du" A d^^’^x*, (2) 

where d”*x = dx^ A ... A dx”^ and d”^“^Xi = i(9/(9x®)d™'x. It is clear from this coordinate 
expression that 12 is a multisymplectic form on Af{T*E). 

As Af‘{T*E) is, in fact, a vector bundle over E, we can consider its quotient by any vec¬ 
tor subbundle. The reduced dual jet bundle of tt, denoted is the quotient of the ex¬ 

tended dual jet bundle, Af'(T*E), by constant affine transformations along the fibers of tt^, 
and is diffeomorphic to the quotient bundle Af" (T*E)/Af'(T*E). The natural quotient map is 
p: Af^(T*E) J^TT*. 

It can be proved that J^tt* may be endowed with the structure of a smooth manifold and, 
moreover, p: Af(T*E) J^tt* is a smooth vector bundle of rank 1. In addition, we have the 
canonical projections : J^tt* —>■ E and vf^: J^vr* —>■ M. 
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Finally, adapted coordinates in E induce coordinates in such that 

the coordinate expression of the natural quotient map is 

where are the induced coordinates in A™(T*£'). In these coordinates, the natural 

projections are given by 

= ; 7r^(x*, = (x*). 

2.3 Multivector fields 

In this section we give a short review on multivector fields and their relation with integrable 
distributions (see [25] for details). 


Locally decomposable multivector fields. Integrability conditions 

A multivector held of degree k, or k-multivector held, on a m-dimensional smooth manifold M 
is a section of the bundle A^(TM) —)• M, that is, a skew-symmetric contravariant tensor. The 
set of all multivector fields of degree /c in M is denoted A^(M). 

In general, given a /c-multivector held X G for every p ^ M there exists an open 

neighborhood Up C M and Xi,... ,Xr € X{Up) such that 

A’= ^ A...AX,,, 

with G C°°([/p) and k ^ r ^ dimM. If for every p we have r = k, that is, there exists an 

open neighborhood Up U M and Xi,..., X^ G X{Up) such that A” = Xi A ... A on Up, then 
we say that the multivector held X is locally decomposable. 

Let P be a A:-dimensional distribution in M, that is, a /c-dimensional subbundle of TM. It 
is clear that sections of M are A:-multivector helds in M, and that the existence of a 

non-vanishing global section of M is equivalent to the orientability of the distribution V. 

Then, we say that a non-vanishing multivector held X G A^(M) and a A:-dimensional distribution 
T) gTM are locally associated if there exists a connected open set U U M such that X\jj is a 
section of A^DI^. 

As a consequence of this we can introduce an equivalence relation on the set of non-vanishing 
/c-multivector helds in M as follows: two /c-multivector helds X,X' G X^{M) are related if, and 
only if, they are both locally associated, on the same connected open set U C M, with the same 
distribution V. In addition, in this case there exists a non-vanishing function / G C°°{U) such 
that X' = fX on U. The equivalence classes of this quotient set will be denoted by {X}ij. 
Then, one can prove that there is a bijective correspondence between the set of /c-dimensional 
orientable distributions "D C TM and set of equivalence classes {X}m of non-vanishing, locally 
decomposable /c-multivector helds in M. 

If A” G X^{M) is a non-vanishing, locally decomposable /c-multivector held and 1/ C M is 
a connected open set, then the distribution associated to the equivalence class {X}if will be 
denoted by Vu{X). If 17 = M, then we write simply 'D{X). 

A non-vanishing, locally decomposable multivector held X G X^{M) is said to be integrable 
(resp. involutive) if its associated distribution 'Du{X) is integrable (resp. involutive). It is clear 
then that if A” G X^{M) is integrable (resp. involutive), then so is every other in its equivalence 
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class {Al}, and all of them have the same integral manifolds. Moreover, Frobenius theorem 
allows us to state that a non-vanishing and locally decomposable multivector held is integrable 
if, and only if, it is involutive. Nevertheless, in many applications we have locally decomposable 
multivector helds fb € X^{M) which are not integrable in M, but integrable in a submanifold 
of M. A (local) algorithm for Ending this submanifold has been developed [25] . 


Multivector fields in fiber bundles and jet bnndles. Holonomy condition 

We are interested in the particular situation of a fiber bundle and, more precisely, of jet bundles. 

Let TT: Li —M be a fiber bundle, with dim M = m and dim E = m n. A multivector field 
X € is said to be n-transverse if at every point u £ E we have {i{X){7r*Lo))y ^ 0 for 

every uj € 17™(M) satisfying uj{'K{y)) ^ 0. It can be proved that if X is integrable, then the 
TT-transverse condition is equivalent to requiring the integral manifolds of X to be local sections 
of vr. In this case, \i (j): U F M ^ E is a local section with = u and is the integral 
manifold of X, then T„(Im(/>) = 'Du{X). 

Now, let us consider the first-order jet bundle of vr, J^vr. A multivector field X £ X'^{J^'k) 
is holonomic if X is integrable, vr^-transverse, and the integral sections of X are holonomic. 

As in mechanics, the holonomy of a multivector field may be characterized using the geometry 
of J^TT. First, a vr^-transverse and locally decomposable multivector field X £ X™'(J^7r) is said 
to be semi-holonomic, or a SOPDE multivector field if, and only if, i{0)X = 0, where 0 is 
the canonical structure form in J^vr. Then, the relation between integrable, holonomic and 
semi-holonomic multivector fields in J^vr is given by the following result from [25] . 

Theorem 1. A multivector field X £ X™'(J^7r) is holonomic if, and only if, it is integrable and 
semi-holonomic. 


In natural coordinates, let X £ X™'(J^7r) be a locally decomposable and vr^-transverse mul¬ 
tivector field. From the results in [25] . this multivector field X may be chosen to have the 
following coordinate expression 



with / being a non-vanishing local function. Then, the condition for X to be semi-holonomic 
gives the mn equations /" = u^. In addition, from the results in [25], we know that the 
necessary and sufficient condition for a locally decomposable multivector field to be integrable 
is that its associated distribution is involutive, which is equivalent to requiring the m(m — l)/2 
conditions \Xj,Xk] = 0, with 1 ^ j < k ^ m. In coordinates, these gives the following system 
of nm{m? — l)/2 partial differential equations for the component functions F}" 


iW, - Ffi, = 0 


dF, 


k,j 


k,i 


dxl 


-|- n' 


dud 




dF°^- 

dx^ 


,dF: 


- w. 




dud 


-F, 


dFl 




kl 




= 0 . 


(3) 


Remark. It is important to point out that a locally decomposable, vf^-transverse and semi- 
holonomic multivector field X may not be holonomic, since the SOPDE condition is not a 
sufficient nor necessary condition for the multivector field to be integrable. On the other hand, 
the integrability of a multivector field does not imply that the integral sections are holonomic: 
as in mechanics, a multivector field may admit integral sections through every point in J^tt, but 
these integral sections may not be projectable to the base manifold. 0 
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Relation with jet fields 

Let vr: E ^ M he a fiber bundle, with dim M = m and dim E = m n, and J^vr the first-order 
jet bundle of vr. 

Definition 6. A jet field in E is a global section of the projection tt^ : — >■ E. 

It is proved in m that there is a bijective correspondence between jet fields in E and 

connections V € A\{T*E) ( 8 ) X{E). Therefore, there is a bijective correspondence between jet 
fields in E and distributions in E. We denote the unique distribution in E associated to 

the jet field T G r( 7 r^). This enables us to give the following definitions. 

Definition 7. Let T G r(7r^) be a jet field and its associated distribution in E. 

1. T is said to be orientable if, and only if, T>('I') is an orientable distribution in E. In 

particular, if M is orientable, then every jet field is also orientable. 

2. T is said to be integrable if, and only if, D(T) is an integrable distribution. 

3. A section 4> G r(7r) is an integral section of T if, and only if, o = j^fi. In particular, 
T is integrable if, and only if, it admits integral sections through every point of E. 

With these notations, the relation between multivector fields and jet fields is given by the 
following result, stated in [25] . 

Theorem 2. There is a bijective correspondence between the set of orientable jet fields T G r(7r^) 
and the set of equivalence classes of locally decomposable and ir-transverse multivector fields 
{T} C X^{E). They are characterized by the fact thatT>{J>) = T>{X). In addition, the orientable 
jet field T is integrable if, and only if, so is every X in the equivalence class. 


3 The Hamilton-Jacobi problem in the Lagrangian formalism 


The geometrical setting for the Lagrangian formalism for multisymplectic field theories is the 
following (see, for instance, [21 EH EH EH |30l Ell ESI [53] for more details). Let vr: S ^ M be a 
fiber bundle modeling the configuration space, where M is a m-dimensional orientable smooth 
manifold with fixed volume form r] G and dimS = m -|- n. Let L G f 7 ™'(J^ 7 r) be a 

Lagrangian density containing the physical information of the theory, that is, a vf^-semibasic 
m-form. We denote L G C°°(J^ 7 r) the function satisfying L = L{Tt^)*rj, which we call the 
Lagrangian function associated to E and p. Using the canonical vertical endomorphism v G 
r(T*J^ 7 r ( 8 )ji 7 r TM ( 8 )ji,r ^(^^))! the Cartan forms Qc = i{r')dE + T G and Qc = 

—d0£ G are constructed, with coordinate expressions 


r)T 

Qc = A d^-^Xi - 

du? 


u°‘ - L \ d'^x 
duf ^ X, 


a2 T o2 T 

^c = ^ . du'^ A du^ A d^-^Xi + —-^ du“ A du^. A d^-^Xi 

dufduP duTduj ^ 


+ \K 


(92 L 


dufduh dud 


^ du^ A d”^x + dw^ A d™x . 


dufduj 


(4) 

(5) 


Then, the Lagrangian problem for first-order classical field theories is the following: to find 
a m-dimensional, -transverse and integrable distribution T>c in J^vr such that the integral 
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sections ofVjr are prolongations of sections (j) € r(7r) satisfying 

i{^)^C = 0 ) for every X E jC( J^vr). (6) 

If the Lagrangian density is regular, then the Cartan (m + l)-form ilc is multisymplectic, and 
then there exists such a distribution, although it is not necessarily integrable. In the following 
we assume that the Lagrangian density L is regular, and that the distribution TJn is, in addition, 
integrable. 

From the results in Section 12.31 this distribution TJn is associated with a class of holonomic 
multivector fields {Xc} C X™'(J^ 7 r) satisfying the equation 

i{Xc)^c = 0, for every Xc E {Xc} ■ (7) 

The same comments apply in the regular case: if the Lagrangian density is regular, then there ex¬ 
ists a class of multivector fields {Xc\ C X^(J^7r) solution to equation ([7|) which is vr^-transverse 
and SOPDE, but not necessarily integrable. In the following we assume that the Lagrangian 
density is regular and that every multivector field in the class is, in addition, integrable. This 
class is denoted by {T/;} along this work. 


3.1 The generalized Lagrangian Hamilton-Jacobi problem 

Following the patterns in m, we first state a generalized version of the Hamilton-Jacobi problem 
in the Lagrangian formalism. 

Definition 8 . The generalized Lagrangian Hamilton-Jacobi problem consists in finding a jet 
field 'L E F( 7 r^) and a m-dimensional and integrable distribution H in E such that if y G r( 7 r) 
is an integral section ofT>, then 'L o 7 E F( 7 r^) is an integral section ofT>c, that is, 

Tulm{y)=Vu Vu E Im(7) TsIm(T o 7) = V u E Im(^'o 7). (8) 


From the results in Section 12.31 since both T> and T>£ are associated with their corresponding 
classes of multivector fields, the problem can be stated equivalently in terms of multivector 
fields as the search of a jet field T E F( 7 r^) and a class of locally decomposable and integrable 
multivector fields {X} C di'^{E) such that if 7 E F( 7 r) is an integral section of every multivector 
field X E {X} then o 7 E F( 7 f^) is an integral section of every multivector field Xc E {Xc} 
solution to equation ([7]), that is, 

Xoy = A ™7 VT E {X} ^ Xco{Jfoy)= ^Xc E {Xc} , (9) 

where M —>■ A™(TFi) denotes the canonical lift of 7 to A”*(T£'). In the following we 

denote by {X} C X^{E) the class of locally decomposable and integrable multivector fields 
associated with the integable distribution T> m. E. The diagram illustrating this equivalent 
formulation of the generalized Lagrangian Hamilton-Jacobi problem is the following 


A’"(^o7) 



A’"7 


7 


A" 


A”*(tf;) 
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where the interpretation is: if the lower diagram formed by 7, X and is commutative for 
every X € {X}, then the upper diagram formed by '£' 07 , Xc and 07 ) is also commutative 

for every Xc G {Xc}- 

Remark. Since the section 'L o 7 G r(7r^) is an integral section of TJc (or, equivalently, an 
integral section of the associated class of holonomic multivector fields), in particular it must 
satisfy equation (l6|), that is, 

('h o 7 )* i{X)Q.c = 0 , for every X € j£( J^vr). 

Nevertheless, observe that the action of the m-form o j)* i[X)Qc G on m tangent 

vectors Vi G Tj,M, with x G M, is defined as 

((^ o 7 )* i{X)nc), {vu...,Vm) = (L!£)v,(-,(.))(X(d/( 7 (x))), T,(d/ o 7 )(ui), ..., T,(d/ o j){vm)) , 

from where we observe that X('I'( 7 (x))) G Im('I'o 7 ) c Im('I'), that is, the vector 

field X is tangent to the submanifold Im(^') ^ J^vr. Therefore, in this particular situation, 
equation ([ 6 ]) is equivalent to 

(T o 7 )* i{X)Qc = 0 , for every X G A( J^vr) tangent to Im(T). (10) 

0 

Remark. Since every integral section of the distribution He is the prolongation of a section of 
TT, this holds, in particular, for the section Toy, and we have T o 7 = j^cj) for some (j) G r( 7 r). 
Now, composing this last equality with the natural projection vr^, we obtain y = (f. Then, 
replacing by 7 in the previous expression, we have T o 7 = from where we deduce that 
if H is an integrable distribution, then the jet field T G r( 7 r^) is integrable, and every integral 
section of H is an integral section of T. Moreover, this enables us to reformulate the generalized 
Lagrangian Hamilton-Jacobi problem as follows: 

The generalized Lagrangian Hamilton-Jacobi problem consists in finding an integrable jet field 
T G r(7r^) such that if y ^ r(7r) is an integral section 0/T, then j^y G r(7r^) is an integral 
section ofHc- 

Nevertheless, we stick to the statement in Dehnition [ 8 l or the equivalent formulation given in 
terms of multivector helds, in order to give several equivalent conditions to being a solution to 
the generalized Lagrangian Hamilton-Jacobi problem. 0 

It is clear from this last remark that the distribution H in E, the jet held T G r(7r^) and the 
distribution He in J^tt are closely related. In fact, we have the following result. 

Proposition 2. The jet field T G r(7r^) and the distribution H in E satisfy condition ([8]) if, 
and only if, H and He are -related, that is, for every Xe G {Xe] (resp., for every X G {X}) 
there exists X G {X} (resp., Xe G {Xe}) such that Xe o T = A^TT o X. 

Proof. We prove this result in terms of the associated classes of multivector helds. Let 7 G r(7r) 
be an integral section of H, which is equivalent to 7 being an integral section of every X G {X}, 
and let Xe G {Xe} be a representative of the equivalence class. Then we have 

Xeoj!oy = A™(¥t^) = A"*TT o = A™TT o A” o 7 , 

where the multivector held X G {X} in the last equality exists since 7 is an integral section of H. 
Then, since X is integrable, it admits integral sections through every point in E, and therefore 
we have proved that for every Xe G {Xe} there exists X G {X} such that Xe o T = A”^T'I' o X. 
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Reversing this reasoning we prove that for every X € {X} there exists G {Xc} such that 
o 'L = o X. Therefore, the distributions V and Vc are T-related. 

Conversely, let us suppose that T> and T>c are T-related. Then, if 7 € r( 7 r) is an integral 
section of T) and Xc € {Xc} is a representative of the equivalence class, we have 

Xcoj!oy = A’^TT 0^07 = A™T^ o A ™7 = , 

where X G {X} in the first equality is any multivector field in {X} which is T-related to the 
given Xc- That is, the jet field T and the distribution V satisfy condition ([ 8 ]). ■ 

A straightforward consequence of this last result is the following. 

Corollary 1. If the jet field T G r(7r^) and the distribution V in E satisfy condition ([ 8 ]), then 
for every X G {X} there exists a multivector field Xc G {Xc} such that X is given by 

A” = A^^Ttt^ oXco^ . 

Proof. Let X G {X} be an arbitrary representative of the equivalence class. From Proposition [2] 
we know that if T and T> satisfy condition ([81), then there exists a multivector field Xc G {Xcf} 
which is 'L-related to the given A, that is, Xc o T = A^TT o X. Then, composing both sides 
of this last equality with the map A^Ttf^: A”^(TJ^ 7 r) A"*(T£'), and bearing in mind that 

T G r( 7 r), we obtain 

A^^Tvri oXcoJf = A^Ttt^ o A™TT o A = A"*T( 7 r^ o T) o A = A . 


That is, every multivector field A G {A} is completely determinated by the jet field T G 
r( 7 r^) and some multivector field Xc G {Xc}-, and it is called the multivector held associated to 
T and Xc- The diagram illustrating this situation is the following: 


J^TT-^^ 



E -^ A^{TE) 


In particular, the distribution D is completely determinated by the jet field T G r(7r^) and the 
distribution Vc as for every u G Im('i;') --G J^vr, or V = T7r^(T>/:|jj^(,^)), 

and it is called the distribution associated to T. 

Remark. From Corollary [1] we deduce that if the jet field 'L is integrable, then so is every 
multivector field A G {A}, and hence the distribution V. Moreover, taking into account the 
remark in page dH we deduce that the jet field T is integrable if, and only if, the distribution 
P is integrable, and they have the same integral sections. Hence, T and P are associated in the 
sense of [25] fSection 12.31) . that is, they dehne the same horizontal subbundle of TE. 0 

Taking into account Corollary [T] it is clear that the search for a jet field F( 7 r^) and a distri¬ 
bution P in P satisfying condition ([ 8 |) is equivalent to the search of a jet field T satisfying the 
same condition with the associated distribution {V C TP. Therefore, we can give 
the following definition. 
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Definition 9. A solution to the generalized Lagrangian Hamilton-Jacobi problem is an inte- 
grable jet field 'L € r( 7 r^) such that if y ^ r( 7 r) is an integral section of the m-dimensional and 
integrable distribution in E, then 'I' o 7 G r( 7 f^) is an integral seetion ofVc- 

Now we state the following characterizations for a jet field to be a solution to the generalized 
Lagrangian Hamilton-Jacobi problem. 

Proposition 3. Let 'L G r( 7 r^) he an integrable jet field. Then, the following statements are 
equivalent. 

1. ^ is a solution to the generalized Lagrangian Hamilton-Jacobi problem. 

2. The distribution He in J^tt is tangent to the submanifold Im('I') ^ J^vr, that is, {Hc)u ^ 
TuIm('L) for every u G Im('I'). 

3. 'L satisfies the equation 


7* = 0 , for every Y G ?C{E ), 

where 7 G r(7r) is an integral seetion of the associated distribution 


Proof. 

(1 2) Assume that 'h G r( 7 r^) is a solution to the generalized Lagrangian Hamilton- 

Jacobi problem, and let 7 G r( 7 r) be an integral section of the integrable distribution H = 
Then, since T and H satisfy condition ([ 8 |) and it is clear that Im('I' o 7 ) C 
Tu Im(T) for every u G Im(Jt 07 ), we have 

= Ts Im(T o 7 ) C Ts Im(T) V tZ G lm(T o 7 ). 

Finally, since T is integrable, for every u G Im('I') there exists an integral section 7 G r( 7 r) such 
that u G Im(T o 7 ), and therefore we have proved {Hc)u C Ti 2 lm('l/) for every u G Im('I'). 

For the converse, assume that the distribution He in J^vr is tangent to the submanifold 
Im('L) J^TT, that is, we have {He)u Y T^Im('I') for every u G Im('I'), which is equivalent to 
(He)q,(u) ^ Im('I') for every u € E. We deduce from this that for every w G there 

exists a Vu & T^E such that w = T„T(uu). Hence, for every u €z E we define a m-dimensional 
subspace H^ C T^E as follows 

Hu = {v e TuE I w = T„T(u) , w G {He)^{u)} ■ 

Then, we define the m-dimensional distribution H va. E as H = 
a smooth and integrable distribution, since it satisfies H = T 7 r^(P£|jj^^^j), 
are smooth and integrable. Now, let 7 G F( 7 r) be an integral section of H. 
of H, the condition for 7 to be an integral section of H gives 

T„Im( 7 ) = T,j,(„) 7 r^((T>£),i,(„)) \/u G fm( 7 ). 

Composing this equality with the map T^T: T„Fi —)> T,j,(„)J^ 7 r, and bearing in mind that 
T o TT^ = Idjj^(^), we have 

{T^ch{u) = T„T(T„fm( 7 )) = Tq>(„) Im(T o 7 ) Vu G Im( 7 ), 
which is clearly equivalent to 

= Ti 2 lm(T o 7 ) ViiGfm(To 7 ), 


It is clear that H is 
and both He and T 
Then, by definition 
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that is, to condition ([S]). Thus, 4^ is a solution to the generalized Lagrangian Hamilton-Jacobi 
problem. 

(1 3) Let 41 € r( 7 r^) be a jet field solution to the generalized Lagrangian Hamilton-Jacobi 

problem, and 7 € r(7r) an integral section of Then the section 4 ^ o 7 G r(7r^) is 

an integral section of In particular, T o 7 is a solution to equation ([U]), that is, 

( 4 / o 7)* i(X)Qc = 0 , for every X € T( J^tt) . 


Calculating, we have 


(T o 7 )* i(X)nc = i(X)nc) = 7 * , 

where V G X(-E) is a vector field 4'-related with X. Nevertheless, since equation ([ 6 |) holds for 
every vector field in J^tt, we have proved 

7 * z(y)(4'*H£) = 0 , for every Y € X(E). 

For the converse, let 7 G r( 7 r) be an integral section of the distribution T 7 r^(T>£|jj^(^)), and 
T G r( 7 r^) an integrable jet field. By the hypothesis we have that 7 * i(y)(4'*H£) = 0 for every 
y G X{E). Then, computing, we have 

7 * i(y)( 4 /*L!£) = 7 ^( 4 /* iix)nc = ('h o 7 )* i{x)nc, 

where X G X(J^ 7 r) is a vector held T-related to Y. In particular, X is tangent to Im(4') ^ J^vr 
by construction, and we have proved that T o 7 G r( 7 r^) is a solution to equation 

(4' o 7 )* i{X)Q.c = 0 , for every X G X(J^ 7 r) tangent to Im(4'), 

that is, to equation (fTOjl . which is equivalent to equation Q for a section of the form Toy. 
Therefore, the section T o 7 G r( 7 r^) is an integral section of T>c-, and therefore 4* is a solution 
to the generalized Lagrangian Hamilton-Jacobi problem. ■ 


A straightforward consequence of the above result is the following. 

Corollary 2. Let 4* G r(7r^) he an integrable jet field solution to the generalized Lagrangian 
Hamilton-Jacobi problem. Then the integral sections ofT>c with boundary conditions in Im(4') 
project to the integral sections of'D = T 7 r^(T>£|jj^(^)). 

Proof. Let 4* G r( 7 r^) be a solution to the generalized Lagrangian Hamilton-Jacobi problem, 
and ifc £ r( 7 r^) an integral section of He with boundary conditions in Im(4'). Then, since the 
distribution He is tangent to Im(4') by Proposition [21 we have that ^ Im(4'), and hence 

T„Im( 7 r^ o V’£) = Tv,>(„) 7 r^(Tq>(„) Im(V’/:)) = Tvi/(„) 7 r),;,(„)) = , 

where we have used that Im(V’£) ^ Im(4') and 4' o vr^ = Id|jjjj(-,j,^. ■ 


Coordinate expression. Let (x*), 1 ^ i ^ m, be local coordinates in M such that rj = d^'x = 
dx^ A... A dx™', and (x*, u"), 1 ^ a ^ n local coordinates in E adapted to the bundle structure. 
Then, the induced coordinates in J^vr are (x*,tt",u“), which coincide with the local coordinates 
adapted to the bundle structure tt^ : J^vr — )■ E. In these coordinates, a jet field 4' G r(7r^) is 
given locally by 4'(x*,ri") = (x*, rt", V'f), where (x*, rt“) are local smooth functions on E. 
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Let us compute the local condition for a jet field 'L € r( 7 r^) to be a solution to the generalized 
Lagrangian Hamilton-Jacobi problem. From Proposition [3] we know that this is equivalent to 
require the distribution Vc in J^vr to be tangent to the submanifold Im(J;') J^vr, or, in terms 
of the class of multivector fields {<L/:} C X™'(J^ 7 r) associated to to require every multivector 
field in the class to be tangent to Im('I'). From [25] we know that a representative G {<L£} 
which is locally decomposable, vr^-transverse and semi-holonomic may be chosen to have the 
following coordinate expression 



where the functions are the solutions to the Euler-Lagrange equations 

dL d‘^L g d‘^L a d‘^L , , 

du°‘ dufdx'^ * dufdul^ ’ 

^ J 

in addition to the integrability conditions (|3|) (if necessary). Observe that every multivector field 
in the class {Xc} is obtained by multiplying this representative by an arbitrary non-vanishing 
function / G C°°(J^ 7 r). Then, bearing in mind that the submanifold Im('I') ^ J^vr is locally 
defined by the mn constraints V'f “ = 0 ) the condition for this particular X^ to be tangent 

to Im('I') gives the following partial differential equations 


dxl ^ du‘^ 



= 0 . 

Im('I') 


( 12 ) 


This is a system of nm? partial differential equations with nm unknown functions '0^, that is, 
we have more equations than unknown functions. 

Remark. Recall that the n equations (jlip do not enable us to determinate all the m^n coefficient 
functions and, in general, there are — arbitrary functions. Therefore, equations (1121) 

may fix not only the coefficients V’f of the jet field ik, but also some of the remaining functions 
F^j of the Euler-Lagrange multivector fields which are solutions to the field equation ([7|). In 
this way, we have a system of m'^n partial differential equations with n{m?‘ + m — 1) unknown 
functions. Note that, even in the most favorable cases, there still are n(m —1) arbitrary functions 
to be determined, which may be fixed by the integrability condition ([3|) or not. 0 


Remark. On time-dependent mechanics, that is, for m = 1 we obtain exactly n partial differ¬ 
ential equations and n unknown functions, since there are no arbitrary functions on the vector 
field solution to the Lagrangian dynamical equation. 0 


3.2 The Lagrangian Hamilton-Jacobi problem 

As in mechanics (see mm), to solve the generalized Lagrangian Hamilton-Jacobi problem is, 
in general, a very difficult task, since it amounts to find {m + n)-dimensional submanifolds of 
J^vr such that the m-dimensional distribution T>c is tangent to them. Because of this, we impose 
an additional condition on the jet field 'k G F(7r^) in order to consider a less general problem. 

Definition 10. The Lagrangian Hamilton-Jacobi problem consists in finding a jet field T G 
F( 7 r^) solution to the generalized Lagrangian Hamilton-Jacobi problem satisfying that = 0. 

Such a jet field is called a solution to the Lagrangian Hamilton-Jacobi problem. 

With this new assumption we can state the following result, which is a straightforward 
consequence of Proposition [3l Corollary [2] and the results in [ 12 ]. 
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Proposition 4. Let 4' € r(7r^) be an integrable jet field satisfying = 0. Then the following 

statements are equivalent: 


1. ^ is a solution to the Lagrangian Hamilton-Jacobi problem. 

2. The submanifold Im(4') ^ J^vr is m-Lagrangian and the distribution Vc is tangent to it. 

3. The integral sections of Vc with boundary conditions in Im(4') project onto the integral 
sections of V = T7r^{'Dc\i^(^q,)). 

Coordinate expression. In coordinates, we have 
q 2 i 02 r 


= 


Ld^-^Xi 

dufdul^ dufdui du^ ) 


+ 


+ 


"2 '^k 

d‘^L dfil a d'^L 

- - + 'ipi - 




du“ A d™x. 


^dufdu^ * du^du°‘ 9u°‘ du^du^du'^ ^ 

Hence, the condition 'I'*H£ = 0 gives the following system of n(l + m{n — 1)) partial differential 
equations 


d^L 


+ 


d'^L dfii 


dufdu/^ dufdu^. duh 


= 0 


d'^L dfi^ 

dufdui 


+ 


/3 


d^L 


* duldu°‘ 


s dfil 

du^Hui du^ 


= 0 , 


'' “'k ... V. ...J .. 

where l^a<;0^nin the first set, and 1 ^ a ^ n in the second. These two sets 

of equations may be combined to obtain the following system 


d’^L 


+ 


d'^L dfii 


dufdud dufduf. du^ 


= 0 


d'^L dfil dL 


dufdui 


du^ 


= 0 . 


(13) 


On the other hand, observe that since 4'*0£ = —d('I'*0£), the condition 4'*0£ = 0 is 
equivalent to requiring the m-form 4/*0£ G 17"*(Li) to be closed. In particular, using Poincare’s 
Lemma, the m-form 4'*0£ is locally exact, that is, there exists a (m — l)-form uj € 17”*“^([/), 
with U ^ E an open set, such that dtu = 'I'*0£. Moreover, since 0£ is 7r^-semibasic, so is 'I'*0£, 
and therefore u must be 7r-semibasic. In coordinates, bearing in mind the coordinate expression 
of the Cartan m-form, we obtain 

dT / dr \ 

= ^du- A d-^x. - - Lj d-x . 

In addition, the coordinate expression for a generic vr-semibasic local (m — l)-form a; in i? is 

u = W^d^-^Xi , 


where VL* € C°°{E) are local functions. Prom this we deduce the local expression of the m-form 
dw, which is 

™ dW^ dW^ 1 

dw = y y- d"*x + y- dn" A d^-^Xi. 

jL^ Pi.,,a 


i=l 


Finally, requiring do; = 'I'*0£, we obtain 


^dw^ dL 


i=l 


9x* 


* du? 


Im(^) 


- L(x\u“,V’f) = 0 ; 


dW^ dL 


du°‘ duf 


Im(^) 


which may be combined to give the equation 

m 


„ dW^ dW^ 

Z- L{x\u^,fif) = 0, 


2=1 


which is the Hamilton-Jacobi equation in the Lagrangian formalism. 


(14) 
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3.3 Complete solutions 

In the above Sections we stated the Hamilton-Jacobi problem in the Lagrangian formalism, and a 
jet field 'h E r(7r^) solution to this problem gives a particular solution to the Lagrangian problem 
in the form of a submanifold the phase space J^tt. Nevertheless, this is not a complete solution 
to the Lagrangian problem, since only the integral sections of the distribution Vc with boundary 
conditions in Im('I') can be recovered from the solution to the Hamilton-Jacobi problem. 

Hence, in order to obtain a complete solution to the problem, we need to endow the phase 
space J^TT with a foliation such that every leaf is the image set of a jet field solution to the 
Lagrangian Hamilton-Jacobi problem. The precise definition is: 

Definition 11. A complete solution to the Lagrangian Hamilton-Jacobi problem is a local 
diffeomorphism U x E ^ J^ir, with U C an open set, such that for every A E C, the 
map 'La(*) = ‘h(A, •): E —>■ J^vr is a jet field in E solution to the Lagrangian Hamilton-Jacobi 
problem. 

Remark. An alternative, but equivalent, definition of a complete solution consists in giving the 
full set of jet fields {Ta E r(7r^) | A E [/ C M™”} depending on mn parameters, instead of the 
local diffeomorphism <1>. 0 

From the definition we deduce that a complete solution endows the Lagrangian phase space 
J^vr with a foliation transverse to the fibers such that every leaf has dimension m -|- n and the 
distribution He is tangent to it. 

It follows from this last comment that a complete solution to the Lagrangian Hamilton- 
Jacobi problem enables us to recover every integral section of the distribution He solution to 
the Lagrangian problem, that is, we can recover every section solution to the Euler-Lagrange 
equations for classical field theories. In particular, let be a complete solution, and let us 
consider the following set of distributions in E: 

{Pa = C TE I A E 17 C R-"} , 

where 'I'a(*) = ‘h(A,«). Then, the integral sections of H\, for different values of the parameter 
A E 17, provide all the integral sections of the distribution He solution to the Lagrangian problem. 
Indeed, let jl,4> E J^vr be a point, and let us denote u = (j){x) = Then, since is a 

complete solution, there exists Xo € U such that <I>(Ao,u) = 'I'a^(u) = fl-f, and the integral 
sections of Hx^ through u, composed with Ta^, give the integral sections of He through jlcj). 

4 The Hamilton-Jacobi problem in the Hamiltonian formalism 

As in the Lagrangian formalism, the conhguration bundle in the Hamiltonian formulation for 
multisymplectic classical field theories is a bundle n: E ^ M, where M is a m-dimensional 
orientable manifold with fixed volume form rj E 12”^(M), and dimE = m n. Two phase 
spaces are considered in this formulation: the extended multimomentum bundle A™(T*E) and 
the restricted multimomentum bundle J^n* introduced in Section [2.21 

The physical information is given in terms of a Hamiltonian section h E r(/r), which is 
specihed by a local Hamiltonian function H E C°°(J^'7r*), that is, we have = 

—H,pI^). Then, from the canonical forms in Af^{T*E), and using this Hamiltonian sec¬ 
tion, the Hamilton-Cartan forms Qh = h*Q E l7"*(J^7r*) and Llh = = —d0/j E 

are constructed, with coordinate expressions 

Qh=pidu‘^ AcT-^Xi-HdTx ; =-dpj, A du" A + dR A d™x . 
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Then, the Hamiltonian problem for a first-order classical field theory is the following: to 
find m-dimensional, rt'^^-transverse and integrable distribution TJh in J^vr* such that the integral 
sections fih € r(7f^) of T>h are solutions to the field equation 

i{X)Qh = 0 , for every X € X( J^vr*). (15) 

Contrary to the Lagrangian formalism for classical field theories, the (m -|- l)-form H/j € 
-g niultisymplectic regardless of the Hamiltonian section h E r(/x) provided, in the 
same way that it occurs in Classical Mechanics in the Lagrangian and Hamiltonian formalisms. 
Therefore, there exists such a distribution P/j, although it is not necessarily integrable. In the 
following we assume that the distribution is integrable. 

From the results in Section 12.31 this distribution is associated with a class of integrable 
and vf^-transverse multivector fields {Xh} C satisfying 

i{Xh)Q.h = 0, for every Xh E {Xh} . (16) 

Same comments apply: since the (m + l)-form Qh is 1-nondegenerate a ff^-transverse solution 
to equation (11611 does exist, but it may not be integrable. In the following we assume that every 
multivector field in the class {Xh} is integrable. 

(For more details on the Hamiltonian formalism of field theories see, for instance, [131 [26l 
EUIMI ESI [53]). 

4.1 The generalized Hamiltonian Hamilton-Jacobi problem 

Following the patterns in [T3] and in previous Sections, we first state a generalized version of 
the Hamilton-Jacobi problem in the Hamiltonian formalism. 

Definition 12. The generalized Hamiltonian Hamilton-Jacobi problem consists in finding a 
section s E and a m-dimensional and integrable distribution H in E such that if y ^ 

is an integral section ofT), then s 07 E F( 7 fg) is an integral section ofDh, that is, 

T„Im( 7 )=P„ Vtt E Im( 7 ) T[^] Im(s 07 ) = (P;j)[^] V [w] E Im(s o 7 ). (17) 

As in the Lagrangian formulation stated in Section El from the results in Section [ 2 . 31 we know 
that both V and Vh are associated with their corresponding classes of multivector fields. Thus, 
we can state the generalized Hamiltonian Hamilton-Jacobi problem in an equivalent way in terms 
of multivector fields as the search of a section s E F( 7 rg) and a class of locally decomposable 
and integrable multivector fields {X} C X™'(H) such that if 7 E F( 7 r) is an integral section of 
every multivector field X E {X}, then s o 7 E F( 7 f^) is an integral section of every multivector 
field Xh E {Xh} solution to equation (fT 6 ]l . that is, 

Xoy = A ™7 yX e {X} ^ Xho soy = A'"(s^) yXh E {Xh} . (18) 

Again, as in the Lagrangian formalism we have the following diagram illustrating this equivalent 
formulation of the generalized Hamiltonian Hamilton-Jacobi problem 

‘xh 
J^TT* 

"'(I’'' 

^ E -^^ A”* (TP) 

A’"7 
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where the interpretation is the same as in the corresponding diagram for the Lagrangian for¬ 
malism; if the lower diagram formed by y, X and A ™’7 is commutative for every X G {X}, then 
the upper diagram formed by s o'y, Xh and A™'(s o 7 ) is also commutative for every Xh € {Xh}. 

Remark. Analogously to the Lagrangian formulation, the section s o 7 G r( 7 r^) is an integral 
section of the distribution TJh solution to the Hamiltonian problem and, therefore, it is a solution 
to the equation (fT^ . that is, 

(s o 7 )* i[X)Q,h = 0 , for every X G jC( J^vr*). 

Then, bearing in mind that the action of the m-form (s o y)* i(X)i}h € 17™'(M) on m tangent 
vectors Vi gTxM (x G M) is defined as 

{{s O 7 )* i{X)Qh)^ {vi,...,Vm) = {^h)s{'y{x)) (silix))) ,Tx{s o y){vi), Tx{s O y){vm )), 

from where we observe that Ai(s( 7 (x))) G Im(s o 7 ) c Im(s), that is, the vector 

field X is tangent to the submanifold Im(s) ^ Hence, we conclude that in this case 

equation (fT^ is equivalent to 

(s o 7 )* i(X)Qh = 0 , for every X G A( J^vr*) tangent to Im(s). 


0 

It is clear from the Dehnition that the distribution H in E, the section s G T{TT'ff) and the 
Hamiltonian distribution TJh in J^tt* are closely related. In fact, we have the following result, 
which is the analogous to Proposition [2] in the Hamiltonian formalism. 

Proposition 5. The section s G r( 7 rg) and the distribution V in E satisfy condition (I17h if, 
and only if, T> and 1)^ are s-related. 

Proof. This proof follows exactly the same patterns as the proof of Proposition [5j ■ 


As in the Lagrangian formalism, a straightforward consequence of this last result is the 
following. 

Corollary 3. If the section s G r(vrg) and the distribution V in E satisfy condition (fT71) . then 
for every X G {X} there exists a multivector field X^ G {X^} such that X is given by 

A” = A”^T7r^ oXhOs. 


Proof. The proof follows the same patterns as in the proof of Corollary [TJ 


That is, the class of integrable multivector fields {X} C XX[E) is completely determinated 
by the section s G r(7rg) and the class of multivector helds {Xh] ^ solution to the 

equation ()16p . and every multivector field X = A^Tyr^ o XhO s a {X} is called the multivector 
held associated to s and Xh. The diagram which illustrates this situation is the following: 


J^TT* -^^ A”*(TJV*) 



E -^ A”^(tl;) 
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As a consequence, the integrable distribution D in i? is completely determinated by the sec¬ 
tion s € r( 7 r^) and the distribution Vh solution to the Hamiltonian problem as L^,r^([a;]) = 
for every [ca] € or V = and it is called the distribution 

associated to s. 

From Corollary [3] we deduce that the search for a section s € r(7rg) and a distribution H in 
E satisfying condition m is equivalent to the search of a section s G r(7rg) such that condition 
dlTD is satisfied with the associated distribution Therefore, we can give the 

following definition. 

Definition 13. A solution to the generalized Hamiltonian Hamilton-Jacobi problem is a section 
s G r(7rg) such that if y ^ r(^) integral section of the m-dimensional and integrable 

distribution then s o 7 G r(7r^) is an integral section ofT>h- 

Proposition 6 . Let s G he a section. Then, the following conditions are equivalent. 

1. s is a solution to the generalized Hamiltonian Hamilton-Jacobi problem. 

2. The distribution in J^vr* is tangent to the submanifold Im(s) ^ that is, (P/i)[a;] C 

T[^] Im(s) for every [w] G Im(s). 

3. s satisfies the equation 


7 * i{Y)d{h o s) = 0 , for every Y G X{E ), 
where 7 G r( 7 r) is an integral section of the associated distribution 

Proof. This proof follows the same patterns as the proof of Proposition [3l bearing in mind the 
properties of the tautological m-form 0 G 17™’(A™(T*H)), that is, we have uj*Q = tv for every 
io G fl^{E). Because of this, we have 

s*VLh = s*{h*n) = {ho s)*n = {ho s)*(-d0) = -d(/i O s)*0 = -d(/i o s) , (19) 

and therefore the equation 

(s o 7)* i[X)Qh = 0 , for every X G X{J^tt*) , 


gives rise to equation 


7 * i{Y)d{h o s) = 0 , for every Y G X{E). 


As in the Lagrangian formalism, a consequence of Proposition [ 6 ] is the following result. 

Corollary 4. Let s G be a section solution to the generalized Hamiltonian Hamilton- 

Jacobi problem. Then the integral sections of V^, with boundary conditions in Im(s) project to 
the integral sections of T> = 

Proof. This proof is analogous to the proof of Corollary [2j ■ 


Coordinate expression. Let (x*), 1 ^ i ^ m, he local coordinates in M such that rj = d™'® = 
dx^ A ... A dx™', and (x*, u“), 1 ^ a ^ n local coordinates in E adapted to the bundle structure. 
Then, the induced coordinates in J^ir* are (x*, which coincide with the local coordinates 

adapted to the bundle structure Tr^: J^vr* —)■ E. In these coordinates, a section s G r(7r^) is 
given locally by s(x*,u“) = (x*, rt“, s),), where s^(x*,u“) are local smooth functions on E. 
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Let us compute the local condition for a section s € r(7rg) to be a solution to the generalized 
Hamiltonian Hamilton-Jacobi problem. From Proposition [6] we know that this is equivalent to 
require the distribution Vh solution to the Hamiltonian problem to be tangent to the submanifold 
Im(s) ^ J^TT*, or, in terms of the associated class of multivector fields {Xh} ^ to 

require every multivector field in the class to be tangent to Im(s). From [53] we know that a 
representative Xh G {Jlh} which is locally decomposable and vr^-transverse may be chosen to 
have the following coordinate expression 


Xk = 


^ / d 

' dp. 


dH d 


du^ 


I _ 


where the functions G\ j are the solutions to the equations 




dH 
du°‘ ’ 


( 20 ) 


in addition to the integrability conditions (if necessary), which in this case give the following 
system of nm{m? — l)/2 partial differential equations for the component functions 

d'^H dH d^H d‘^H d'^H dH d'^H d'^H _ ^ 

dxldp^^ du^dp^^ dp^pdp'^ dx^dpi dp^^ du^dpi dp^^dpL 

dxl dp^p du^ dpp dx^ dpp du^ 


Then, bearing in mind that the submanifold Im(s) ^ J^tt* is locally defined by the mn 
constraints —p^ = 0, the condition for this representative of the class to be tangent to Im(s) 
gives the following partial differential equations 


dH 

dxl dpL 



= 0 . 


Im{s) 


( 22 ) 


This is a system of nm? partial differential equations with nm unknown functions that is, 
we have more equations than unknown functions. 


Remark. Recall that equations (12011 do not determinate uniquely all the coefficient functions 
since there are m^n unknown functions ^ and we have only n equations, which implies 
that, in general, there are n{m? — 1) arbitrary functions. Therefore, equations (|22ll could enable 
us to fix some of the arbitrary functions G)^ ^ of the Hamiltonian multivector fields solution to 
the field equation (|16ll . From this point of view, equations (|22p are a system of m?n partial 
differential equations with n(m^ + m — 1) unknown functions. Note that in the most favorable 
cases there still are n{m — 1) arbitrary functions to be determined, which may be hxed by the 
integrability condition (j2ip or not. 0 


Remark. On time-dependent mechanics, that is, for m = 1 we obtain exactly n partial differ¬ 
ential equations and n unknown functions, since there are no arbitrary functions on the vector 
field solution to the Hamiltonian dynamical equation. 0 


4.2 The Hamiltonian Hamilton-Jacobi problem 

As in mechanics (see |141I18] ). to solve the generalized Hamiltonian Hamilton-Jacobi problem is, 
in general, a very difficult task, since it amounts to find mn-codimensional submanifolds of J^vr* 
such that the m-dimensional distribution is tangent to them. For this reason we require an 
additional condition to the section s G r(7rg), and thus we consider a more particular problem. 
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Definition 14. The Hamiltonian Hamilton-Jacobi problem consists in finding a section s € 
r(7r^) solution to the generalized Hamiltonian Hamilton-Jacobi problem such that s*iJh = 0. 
Such a section is called a solution to the Hamiltonian Hamilton-Jacobi problem. 

Remark. Bearing in mind the properties of the tautological m-form 0 G (A'lfi(T*E)), and 
the calculations in (1191) . the condition s*Qh = 0 is equivalent to the closedness of the m-form 
hos€H^{E). 0 

With this new assumption we can state the following result, which is a straightforward conse¬ 
quence of Propositional Corollary H] and the results on isotropic submanifold of multisymplectic 
manifolds in |12j . 

Proposition 7. Let s G T{7r^) be a section satisfying s*Llh = 0. Then the following statements 
are equivalent: 


1. s is a solution to the Hamiltonian Hamilton-Jacobi problem. 

2. The submanifold Im(s) ^ is m-Lagrangian and the distribution Vh solution to the 

Hamiltonian problem is tangent to it. 

3. The integral sections of TJh with boundary conditions in Im(s) project onto the integral 
sections of T> = T7r^(X>hlim{s))- 


Coordinate expression. In coordinates, we have 


s*VLh = —d(/i o s) = 


dH 

du‘^ 


dHds^p ^dsj 


dpi du°‘ 


2=1 


du" A d"*x + ^ du“ A duH A d”^ 
OuP 


-1, 


Hence, the condition s*Llh = 0 or, equivalently, d{ho s) = 0, gives the following n(l -|-m(n — 1)) 
partial differential equations 


dH dH 54 ^ dsi 

^ p 2=1 


dud du^ ’ 


(23) 


where 1 ^ a ^ n in the first set, and l^Q:</3^nin the second. Therefore, a section s G r(7r^) 
solution to the Hamiltonian Hamilton-Jacobi problem must satisfy the m^n -|- mn{n — 1) -|- n 
partial differential equations ([2^ and ([23l) . 


Remark. Since equations (12311 are the analogous to equations (jl3|l in the Hamiltonian for¬ 
malism we deduce that, in general, equations (j23p may not be locally linear independent 
among themselves or together with (I22p . In particular, we may have less equations than the 
m^n -|- mn(n — 1) -|- n given by (1221) and (f23]l . 0 


Now we recover the classic Hamilton-Jacobi equation for first-order classical field theories. 
Since the Tr^-semibasic m-form h o s is closed, by Poincare’s Lemma it is locally exact, that 
is, there exists a vr-semibasic (m — l)-form io G 17™'“^ (C/), with U C E an open set, such that 
do; = ho s. In coordinates, a semibasic (m — l)-form defined in an open set of E is given by 

CO = W^cT-^Xi , 


with VP* G C°°{E) being local functions. From this, the m-form dw is given by 

5VP* _ 9VP* 


do; = y ^ d"*x + ^ du“ A d-^x, 
dx* 


2=1 


du° 
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Hence, requiring the identity dw = h o s to hold, we obtain 


-Hix\u^,sl) = J2 


i=l 


dW^ 


dW^ 


from where the classic Hamilton-Jacobi equation is deduced 

y ^ + H [ x\ u", 1=0. 


2 = 1 


dx^ 


du° 


(24) 


4.3 Complete solutions 

As in the Lagrangian formalism stated in Section [3l in the previous Sections we stated the 
Hamiltonian formulation of the Hamilton-Jacobi problem for first-order classical held theories. 
As we have proved previously, a section s G r(7r^) solution to the Hamilton-Jacobi problem gives 
rise to a particular set of solutions to the Hamiltonian problem in terms of a submanifold of 
phase space Observe, however, that this section of the bundle J^tt* E is not a complete 

solution to the Hamiltonian problem, since only the integral sections of Vh with boundary 
conditions in Im(s) can be recovered from the solution to the Hamilton-Jacobi problem. 

Therefore, a complete solution to the problem is given in terms of a foliation of the phase 
space J^TT* such that every leaf is the image set of a section solution to the Hamiltonian Hamilton- 
Jacobi problem. 

Definition 15. A complete solution to the Hamiltonian Hamilton-Jacobi problem is a local 
diffeomorphism U x E ^ with U C M”*”' an open set, such that for every A € [/, the 

map sa(*) = E —)• J^vr* is a section of the projection which is a solution to the 

Hamiltonian Hamilton-Jacobi problem. 

As in the Lagrangian formalism, it follows from the dehnition that a complete solution 
provides the manifold J^tt* with a foliation transverse to the fibers satisfying that every leaf is 
(m -|- n)-dimensional and the distribution 1)^ is tangent to it. 

Moreover, from a complete solution to the Hamiltonian Hamilton-Jacobi problem we can 
recover every section solution to the Hamilton-De Bonder-Weyl equations. In particular, let 
be a complete solution, and let us consider the following set of distributions in E 

{Pa = C Ti? I A G [/ C R™-} , 

where sa(») = <h(A,«). Then, the integral sections of Vx, for different values of A G C/, provide 
all the integral sections of the distribution Hh solution to the Hamiltonian problem. Indeed, let 
[w] G J^vr* be a point, and let u = vrg([w]) be its projection to E. Then, since 4> is a complete 
solution, there exists Xq & U such that <h(Ao,u) = s\^{u) = [w], and the integral sections of Vx^ 
through u, lifted to J^vr* by sx„, give the integral sections of through [w]. 

4.4 Equivalence with the Lagrangian formalism 

In previous Sections we have stated the (generalized) Hamilton-Jacobi problem for multisymplec¬ 
tic classical held theories in both the Lagrangian and Hamiltonian formalisms. In the following 
we establish a bijective relation between the solution of the Hamilton-Jacobi problem in both 
formulations in terms of the restricted Legendre map. 

First of all, recall that since the Cartan m-form 0£ G l7”*(J^7r) is vr^-semibasic, we can give 
the following dehnition. 
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Definition 16. The extended Legendre map associated with the Lagrangian C € f7’"(J^7r) is 
the bundle map FC: J^vr —)■ over E defined as follows: 

• • •, • • •, {Tji^Tr^ymrnx))) , 

where Yi € X(J^7r), and therefore Tvr^li G X{E). 


This map verifies tte oFC = vr^, that is, EC is a bundle morphism over E. Furthermore, we 
--—■—'* -- * 

have that EC 0 = 0£ and EC Q = Then, bearing in mind the coordinate expressions ([2]) 
of the tautological m-form 0 G 1?™'(A™(T*£')), and ([H of the Cartan m-form 0£ G 
the coordinate expression of the extended Legendre map is 


EC X* = x~ 


EC = u- 




dL 

duf 


EC p = L-uf 


dL 

duf 


(25) 


The composition of the extended Legendre map EC: A™(T*Fi) with the natural 

quotient map p: Kfif{T*E) —>■ J^vr* gives rise to a bundle morphism poEC: J^vr which 

leads to the following definition. 

Definition 17. The restricted Legendre map assoeiated to the Lagrangian density C G 1?™'( J^vr) 
is the bundle morphism EC: J^tt —>■ J^vr* over E defined as EC = po EC. 

In the natural coordinates of J^vr*, the local expression of the restricted Legendre map is 

dT 

ECE = X* ; ECu^ = ; EC*p^^ = — . 

duf 


As for the extended Legendre map, the map EC: J^vr —>■ J^-ir* satisfies vr^ o EC = vr^, 
EC*E)h = 0£ and EC*Clh = Clc- Moreover, the Lagrangian density C G is regular 

if, and only if, the restricted Legendre map is a local diffeomorphism, and C is said to be 
hyperregular if EC is a global diffeomorphism. Furthermore, in the hyperregular case the 
Hamiltonian section h G F(/i) is defined hy h = ECoEC~^ (or by restriction on the corresponding 
open sets where EC is a diffeomorphism in the regular case). 

Finally, to close the review on the properties of the Legendre maps, we have the following 
fundamental result from [53]. 

Theorem 3. Let (J^7r,£) be a Lagrangian field theory with C G l?™'(J^7r) a hyperregular La¬ 
grangian density, and (J^7r*,H/i) the assoeiated Hamiltonian field theory. 

1. If 4> ^ r(7r) is a solution to equation Q, then the section fih = IFH ° £ r(7r^) is a 

solution to equation m- 

2. Conversely, if fih € r(7r^) is a solution to equation (fTKIl . then the section 4> = ° ^ 

F(7r) is a solution to equation ([6|). 

In particular. Theorem [3] states that the distributions Vc and Of, solutions to the Lagrangian 
and Hamiltonian problems, respectively, are J^£-related and, moreover, it establishes a bijective 
correspondence between the integral sections of both distributions. 

Now we can state the equivalence between the solutions of the (generalized) Lagrangian and 
Hamiltonian Hamilton-Jacobi problems in terms of the restricted Legendre map. First we need 
the following technical results. 
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Lemma 2. Let Ei M and E 2 M he two fiber bundles, F: Ei E 2 a fiber bundle 
morphism, and two F-related k-multivector fields Xi G X^[Ei) and X 2 G 3C^{E2). If si G r(7ri) 
is a section o/vri and we define a section of 7^2 as S 2 = F o si ^ r(7r2), then 

A^^Tvri o Al o Si = A^Tvra o Aa o S 2 G X^{M). 


Proof. As -F: -El —7- -©2 is a fiber bundle morphism (that is, vri = tt 2 o F), and Ai and Aa are 
F-related (that is, A*^T.F o Ai = A 2 o F), we have the following commutative diagram 


A^T-Fi 




A^^T-Fa 


A-i 

El 



X 2 

-Fa 


Then we have 


A’^Ttti o Ai o si = A^T(7r2 o F) o Ai o si = A^'Tvra o A'^TF o Ai o si 
= A^Tvra o Aa o F o si = A^Tvra o Aa o sa • 


Lemma 3. Let Ei Eh,- M and Fa M be two fiber bundles, F: Fi —>■ Fa a fiber bundle 
morphism, and two F-related k-dimensional distributions T>i in Ei and Fa in Fa. If si G r(7ri) 
is a section of tti and we define a section o/vra as sa = F o si G r(7r2), then 

T7ri(Fi|jj^(^^)) = T7r2(F2lim(s2))' 

Proof. Since Fi and Fa are F-related, it suffices to apply Lemma [2] to every pair of F-related 
multivector fields in the classes of multivector fields associated to Fi and Fa. ■ 


Finally, the equivalence Theorem is the following. 

Theorem 4. Let be a Lagrangian field theory with C G i?’"(J^7r) a hyperregular La- 

grangian density, and (J^7r*,F/j) its associated Hamiltonian field theory. 

1. // 'L G r(7r^) is an integrable jet field solution to the (generalized) Lagrangian Hamilton- 
Jacobi problem, then the section s = FC o T G r(7rg) is a solution to the (generalized) 
Hamiltonian Hamilton-Jacobi problem. 

2. Conversely, if s € r(7rg) is a solution to the (generalized) Hamiltonian Hamilton-Jacobi 
problem, then the jet field T = FC~^ o s G r(7r^) is a solution to the (generalized) La¬ 
grangian Hamilton-Jacobi problem. 

Proof. This proof follows the patterns of Theorem 3 in |14j and Theorem 1 in |18] , now in terms 
of distributions. 

In particular, assume that T G r(7r^) is a solution to the generalized Lagrangian Hamilton- 
Jacobi problem. First, if F = T7r^(F£|jj^j-,j,^), F = T7r^(F/i lim(s)) integrable distribu¬ 

tions associated to T and s = FL o T G F(7rg), respectively, then by Lemma[3]we have T> = f>. 
Hence, both distributions are denoted by F. 
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Now we prove that s = J-C o is a solution to the generalized Hamiltonian 

Hamilton-Jacobi problem. Let 7 G r(7r) be an integral section of V. Computing, we have 

7* i(r)d(h os) = j* i(Y)d(h os) = 7* i(V)(s*nk) 

= 7* i{Y){(FC O i{Y)di*{FC*nh) 

= 7* i{Y){j>*nc) ■ 

Hence, since this last expression vanishes for every Y G X{E) by the hypothesis and Proposition 
[3l we have proved 

7* i{Y)d{h o s) = 0 , for every Y G X{E ), 

which, using Proposition O is equivalent to s = TC o vI/ being a solution to the generalized 
Hamiltonian Hamilton-Jacobi problem. 

Now we require, in addition, = 0, that is, '1' is a solution to the Lagrangian Hamilton- 

Jacobi problem. Then we have 

= {FC o = j>*{FC*nh) = d>*nc = o, 

and therefore s = FC o T is a solution to the Hamiltonian Hamilton-Jacobi problem. 

The converse is proved analogously, but using FC~^ instead of FC. ■ 


5 Examples 

5.1 Non-autonomous dynamical systems 

Let us consider the case of non-autonomous dynamical systems, that is, the base manifold M 
is 1-dimensional, and in particular we assume that M = M with the canonical volume form 
rj G 17^ (M). Let us consider a first-order non-autonomous dynamical system with n degrees of 
freedom, and let vr: S —M be the configuration bundle for this system, with dim LI = n-|-l. The 
dynamical information is enclosed on a Lagrangian density C G l?^(J^7r), which is a vr^-semibasic 
1-form. Because of this, let us denote by L G C°°(J^7r) the function satisfying C = L ■ (vr^)*?]. 

Remark. Observe that since M is 1-dimensional, it is either diffeomorphic to the real line M 
or the unit circle S^. The only difference for our calculations is that, contrary to the case of the 
unit circle, M has a canonical global atlas with a distinguished coordinate. Nevertheless, both 
manifolds are orientable and parallelizable, and therefore all the calculations remain the same 
in the unit circle despite the absence of a global atlas. 0 

Local coordinates along this Section are denoted in the usual way for non-autonomous 
dynamical systems. In particular, let {t) denote the global coordinate in R, and let 
1 ^ H ^ n, be local coordinates in E adapted to the bundle structure. Then, the induced 
local coordinates in J^vr, T*E and J^tt* are denoted {t,q"^,v^), it,q^,p,PA) and {t,q^,pA), 
respectively. 


Lagrangian formalism 

The Lagrangian problem for first-order non-autonomous dynamical systems consists in finding 
a class of vf^-transverse and semi-holonomic vector fields {Xc} C X{J^7r) satisfying the field 
equation ([7|), that is, 

i{Xc)Clc = 0, for every Xc G {Xc} ■ 
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If the Lagrangian density is regular, the Cartan 2-form Q G f2^(J^7r) has maximal rank 2n on 
J^vr, and therefore there exists a class of vector fields {X^) solution to the above equation which 
is holonomic. Hence, in the following we assume that the Lagrangian density is regular. 

Remark. Note that the integrability is assured without further assumptions in this setting, 
since every vector field defined on a manifold is always integrable. 0 


The local expressions of the Cartan 1-form 0£ G l7^(J^7r) is 


dL . f dL 


— L] dt, 


(26) 


from where we deduce the coordinate expression of the Cartan 2-form, which is 

,2 r q2j^ 


Vic = 


9 L . Q 

rdgr^ A dq^ + 


dv^dq^ 


+ V 


dv^dv^ 


dq^ A dv^ 


d^L 


dy^dqB Qg 


dL \ n , A d^L 
dq^ Adt + v 


B 


(27) 


B A ^ ^ dv^ A dt. 


dv^dv^ 


Thus, a representative of the equivalence class {Xc} C X(J^7r) of tt^- transverse vector fields 
solution to the above dynamical equation is given in coordinates by 


Xc = f 




A d r^A 9 

. -\- - 1 - - 

\dt dq^ dv^ 


(28) 


where / is a non-vanishing local function, and the n functions are the unique solutions to 
the Euler-Lagrange equations 

dL d^L ^ d‘^L ^ d‘^L 
dq^ dtdv^ ^ dq^dv^ dv^dv^ 


Observe that the vector fields in the class {Xc} VL X(J^7r) are completely determinated (except 
for the non-vanishing function /). In the following we take / = 1 as a representative of the 
equivalence class to state the Lagrangian Hamilton-Jacobi problem. 

For the generalized Lagrangian Hamilton-Jacobi problem, all the Definitions and results 
stated in Section o remain without changes, except for Corollary [2l where “boundary condi¬ 
tions” should be replaced by “initial conditions”. Note that every mention to the integrability 
of the jet field is redundant, since the associated distribution in this case is 1-dimensional, and 
therefore integrable without further assumptions. 

In coordinates, let T G r(7r^) be a jet held locally given by 'I'(t, q^) = (t, q^)), with 

^ being local functions. Then, from Proposition[3]we know that T is a solution to the 

generalized Lagrangian Hamilton-Jacobi problem if, and only if, every vector held in the class 
{Xc} is tangent to the submanifold Im('I') ^ J^tt, which is locally dehned by the constraints 
= 0, 1 ^ H ^ n. Then, the conditions L{Xc){v^ — = 0 give rise to the following 

system of n partial differential equations 


F^\ 

Ilmj-P) 


dip^ 

~dr 






= 0 . 


Now, in order to state the Lagrangian Hamilton-Jacobi problem, we must require in addition 
the jet held 'k G r(7r^) to satisfy the condition = 0. In this situation, DehnitionllOl remains 

unchanged, but the statement of Proposition |4] must be changed as follows. 

















M. de Leon et al.: Hamilton-Jacobi theory in multisymplectic classical held theories 


28 


Proposition 8. Let 'I' € r(7r^) be a jet field satisfying = 0. Then the following statements 

are equivalent: 


1. ^ is a solution to the Lagrangian Hamilton-Jacobi problem. 

2. The submanifold Im('I') ^ J^vr is Lagrangian with respect to the presymplectic form Qc 
and every vector field S {Xc} is tangent to it. 

3. The integral curves of Xc € {Xc} with initial conditions in Im('I') project onto the integral 
curves of X = Ttt^ o Xc o 


Bearing in mind the coordinate expression (I27p of the Cartan 2-form, the pull-back G 

has the following coordinate expression 

d^L . 

^ “ \dv^dqB ^ dv^dv^ dqBJ^^'^ 

f d^L dfi^ d^L dL d^L dfi^ \ ^ 

\dv^dv^ dt dv^dq^ dq^ dv^dv^ dq^) ^ ^ ^ 

Hence, the condition = 0 gives the following system of partial differential equations 

d^L d'^L dfi^ _ _ d^L dfi^ jA 9'^ ^ dL . d‘^L dfi^ _ 

dv^dq^ dv^dv^ dq^ ’ dv^dv^ dt dv^dq^ dq^ dv^dv^ dq^ ^ ’ 

which may be combined to obtain equations (jl3[) in this setting, that is, the following system of 
partial differential equations 

d^L dfi^ d'^L _ d'^L dfi^ dL 

dv^dq^ dq^ dv^dv^ ’ dv^dv^ dt dq^ 


Finally, we state the Hamilton-Jacobi equation in the Lagrangian formalism. Since the 
jet field 'L G r( 7 r^) is a solution to the Lagrangian Hamilton-Jacobi problem, we have that 
d(\J/*0£) = 0. Thus, there exists a local function W G C°°{E) such that dkF = T*©/;. In 
coordinates, the 1 -form 'I'* 0 £ is given by 


9L jy f ^ dL 
dv^ \ dv^ 

Hence, requiring T*©/; = dVF, we obtain 


- L{t,q^,ii^) ) dt. 


dW . 5L ^ . dW dL 


= -r 


+ L{t,q ) ; 


dt dv^ ’ ’ ’ dq^ dv^ ’ 

which may be combined in the following single equation 

dW .dW 

— - Lit,q ,'ip ) = 0 . 


dt 


dq^ 


This is the usual Hamilton-Jacobi equation for a first-order non-autonomous Lagrangian dy¬ 
namical system. 
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Hamiltonian formalism 


Now we state the Hamiltonian formulation of the Hamilton-Jacobi problem in this setting. 
In the natural coordinates of the coordinate expression of the restricted Legendre map 

J-£ : J^TT ^ J^TT* in this setting is the following 


J^£*pa = 


dL 

dv^ 


from where we deduce the coordinate expression of the extended Legendre map, which is 


PA = 


dL 

dv^ 


FC,p = L-v^ 


dL 

dv^ 


Recall that since £ € H^{J^Tr) is assumed to be regular, the restricted Legendre map is a local 
diffeomorphism. For simplicity, in the following we assume that £ is hyperregular, and thus TL is 
a global diffeomorphism (the regular but not hyperregular case is recovered by restriction on the 
open sets where the restricted Legendre map is a diffeomorphism). Let h = TC. o T 
be the Hamiltonian section associated to the Lagrangian density £, and H G C°°(J^7r) a local 
Hamiltonian function which specifies this Hamiltonian section h. 

Using the Hamiltonian section h G r(/i) we define the Hamilton-Cartan forms = h*Q G 
and Llh = h*Ll G f?^(J^7r*), with coordinate expressions 


©h = PAdq'^ — Hdt ; Vlh = dq^ A dpA + dH A dt. 


(29) 


The Hamiltonian problem for first-order non-autonomous dynamical systems consists in find¬ 
ing a class of vf^-transverse vector fields {X^} C X{J^7r*) satisfying equation (fT6|) . that is, 

i{Xh)VLh = 0, for every Xh G {Xh} . 

As in the general setting described at the beginning of Section [H in this formulation the 2-form 
Lth in J^vr* has maximal rank 2n regardless of the Hamiltonian section h G r(^). Therefore, 
there always exists such a class of vector fields. 

Remark. As in the Lagrangian formulation described previously, the integrability is assured 
without further assumptions in this setting. 0 

Then, a representative of the equivalence class {X^} U j£( J^vr*) of vf^-transverse vector fields 
solution to the above dynamical equation is given in coordinates by 

_ J d dH d dH d \ 

^ V dp A dq^ dq^ dp a ) 

As in the Lagrangian formulation, note that the vector fields in the class {X^} C are 

completely determinated (again, except for the non-vanishing function /). In the following we 
take / = 1 as a representative of the equivalence class to state the Hamiltonian Hamilton-Jacobi 
problem. 

For the generalized Hamiltonian Hamilton-Jacobi problem, all the Definitions and results 
stated in Section O remain unchanged, except for Corollary 01 where “boundary conditions” 
should be replaced by “initial conditions”. 

In coordinates, let s G r(7r^) be a section locally given by s{t,q^) = {t,q'^, SA{t,q^)), where 
SA G C°°(£) are local functions. Then, from Proposition [6] we know that s is a solution to the 
generalized Hamiltonian Hamilton-Jacobi problem if, and only if, every vector field in the class 
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{Xh} is tangent to the submanifold Im(s) ^ Since this submanifold is defined locally by 

the constraints pA — sa = 0, 1 ^ A ^ n, we must check if the conditions L{Xh){pA — sa) = 0 
hold, which give rise the following system of n partial differential equations 


dH dsA dH dsA 

dq^ dpB dq^ ’ 


(on Im(s)) . 


Now, in order to state the Hamiltonian Hamilton-Jacobi problem, we must require in addition 
the section s G r(7r^) to satisfy the condition s*ilh = 0. In this case, the statement in Definition 
M and the remark that follows remain unchanged, but the statement of Proposition [7] must be 
changed as follows 

Proposition 9. Let s G r(7r^) be a section satisfying s*Llh = 0. Then the following statements 
are equivalent: 


1. s is a solution to the Hamiltonian Hamilton-Jacobi problem. 

2. The submanifold Im(s) J^tt* is Lagrangian with respect to the presymplectic form Llh 
and every vector field G {Xh} is tangent to it. 

3. The integral curves of Xh G {X^} with initial conditions in Im(s) project onto the integral 
curves of X = Tvr^ o X^ o s. 


In coordinates, bearing in mind the coordinate expression (j29ll of the Hamilton-Cartan 2- 
form D/i, the pull-back s*0,h G H‘^{E) has the following local expression 


= 


dsA dH dH dsB 


+ 


+ 


dt dq^ dpB dq^ 


dg"^ Adt + l^dg^ A dg^ 


Hence, the condition ,s*Qh = 0 gives equations (l23|) . which in this case correspond to the following 
system of partial differential equations 

dsA dH dH dsB dsA dsB 
dt dq^ dpB dq^ ’ dq^ dq^ 


Finally, we state the Hamilton-Jacobi equation in the Hamiltonian formalism. Since the 
section s G T{TT'ff) is a solution to the Lagrangian Hamilton-Jacobi problem, we have that 
/los G Q^{E) is a closed form. Therefore, using Poincare’s Lemma, there exists a local function 
W G C°°{E) such that dIT = ho s. Bearing in mind that the 1-form ho s £ H^{E) is given in 
coordinates by 

ho s = SAdq^ — [H o s)dt, 

the condition ho s = dVF gives the following partial differential equations 


dW 


= -H{t,q'^,SA) 


dW 


dt ’ dq^ 

which may combined to obtain the following single equation 


= SA , 


dW 


+ H 



= 0 . 


This is the classic Hamilton-Jacobi equation for a hrst-order non-autonomous dynamical system 
with n degrees of freedom and Hamiltonian function H. 


















M. de Leon et al.: Hamilton-Jacobi theory in multisymplectic classical held theories 


31 


5.2 Quadratic Lagrangian densities 

Let us consider a classical field theory with n fields depending on m independent variables given 
in terms of a quadratic Lagrangian density, and let vr: i? —>■ M be the configuration bundle, 
with M being a m-dimensional orientable smooth manifold with fixed volume form r] € 1?™'(M), 
and dim Li = m + n. Most of the quadratic Lagrangian field theories can be modeled as follows 
[561 |29]: let us assume that the bundle tt: E —>• M is trivial (that is, E = M x Q, with 
dim <5 = n), so that : J^tt —>■ Li is a vector bundle. Let g he a pseudo-Riemannian metric in 
this vector bundle, L a connection of the projection and V € C°°(E) a potential function. 
Then a quadratic Lagrangian density C € l 7 ^(J^ 7 r) is given by 

In the natural coordinates of J^tt, the m-form L has the following local expression 

d^x, 

where tt") are the coefficients of the metric, which moreover satisfy due 

to the symmetry assumption, and Lf = Lf (x*) are the component functions of the connection. 

In order to simplify the problem, we assume that the connection L is integrable. As a 
consequence, there exist natural charts of coordinates in J^tt such that Lf = 0 , which implies 
that the coordinate expression of the Lagrangian density E reduces to 

/:(x\u",u")= Q 5 >f< + n^*,^“))d-x = L-d-x. (30) 

Observe that the Hessian matrix of the Lagrangian function L associated with L and rj 
coincides with the matrix of the coefficients of g, that is. 



and therefore the Lagrangian density C is regular since 5 is a pseudo-Riemannian metric, and 
in particular non-degenerate. 


Lagrangian formalism 


The local expression of the Cartan m-form Qjr G !?’"(J^vr) obtained from the Lagrangian density 
given locally by (1301) is 


A d—lx, - Q - V{x\ u")) d^x , 

from where we obtain the coordinate expression of the Cartan (m -|- l)-form, which is 


dg, 




A du" A d”^-^Xi - g'-^pdu'^j A du" A d^^-^Xi + g^pU^du^- A d^ 


du^ 


, (dg% p ^ idg^p p dV \ , s . 

\ dx^ ^ 2 du^ ^ du^ ' 


(31) 


( 32 ) 
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Hence, a representative of the class of holonomic multivector fields {Jdc} ^ X"*(J^7r) solution 
to the Lagrangian field equation d?]) is given in coordinates by 



where the functions are solutions to the Euler-Lagrange equations for multivector fields (llip , 
which in this case are 


1 s 

_ _ L _ I! .li -\- _ 

2 ^ 


E 


dxi i 


^9a& P S 7-1/3 ij 


(33) 


in addition to the integrability conditions ([3]) (if necessary). 

First we state the generalized version of the Hamilton-Jacobi problem. Let € r(7r^) be a 
jet field locally given by 'L(x*,u“) = (x*, n“, V'f(x*, n“)), where G C°°{E) are local functions. 
Then, from Proposition [3] we know that T is a solution to the generalized Lagrangian Hamilton- 
Jacobi problem if, and only if, every multivector held Xc € {X^} is tangent to the submanifold 
Im('L) ^ J^TT. Bearing in mind that Im('I') is dehned locally by the constraints = 0, 

the tangency condition gives rise to equations (I12h . where now the component functions are 
solutions to the n equations ([331) . 

In order to state the equations of the Lagrangian Hamilton-Jacobi problem, we must require 
in addition the jet held T G r('?r^) to satisfy the condition = 0. Bearing in mind the 

coordinate expression (j32l) of the Cartan (m -|- l)-form Q,c, the pull-back G is 

given in coordinates by 




'dgi^ 




du^ 


< + 


_ 1 L 

0/3 


du^ A d" A d^ 


m—1 


Xi 


1 f ij ^^j , ij la^^j , ^9ip ,P , 
+ 9sp~^ + 9api’i + 


ld9% 


du^ (9x* ^ 2 du^ 





du-^ A d”^x. 


Hence, the condition = 0 gives the following system of partial differential equations 


8u» ” 


■ ■ dip^ 

Q-P QyS 


= 0 


• ■ dip' 
u _2 
9SP Q^i 


P 


+ 9i,t 


dip 


P 




dg'i 




du^ (9x* ^ 




,/3 »V 

A - ^ 


= 0 . 


which may be combined to obtain equations (|13p . that is, the following system of n(l-|-m(n —1)) 
partial differential equations 


QyS j 0/3 QyS 


p 


= 0 


ij 

95P 


+ 


dg"i 


9x* 


sp I p 

-^j 


1 dfii 


OijS 


2 du^ 




dV 

du^ 


Finally, we deduce the Lagrangian Hamilton-Jacobi equation. Since the jet held T G r(7r^) 
is a solution to the Lagrangian Hamilton-Jacobi problem, we have that T*0£ G Q'^{E) is a 
closed form, and thus there exists a local vr-semibasic (m — l)-form lo G 1?’"“^(F) such that 
do; = In coordinates, bearing in mind the coordinate expression (I3ip of the Cartan 

m-form, we obtain 




©£ = A d 


m—1 


Xi - 


^gip€^i-v{x\un]d^x 
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Then, if the (m — l)-form lo is given locally by a; = ITM™ ^xt, the condition do; = T*0£ gives 
rise to the following equations 


^ CTT* 1 ij a j 

2 ^ ^ + 2 ) = 0 
2=1 


dW^ 

du°‘ 


- W/ 


which may be combined to give the classic Hamilton-Jacobi equation in the Lagrangian formalism 
m, which in this example is 


E 


dW^ 

dx^ 


1 -a0dW^ dW^ 

2 Qy^a QyP 


= 0 . 


(34) 


where denote the coefficients of the inverse matrix of j j which exists since we assume 
g to be non-degenerate. 


Hamiltonian formalism 


Now we state the Hamiltonian formalism for the Hamilton-Jacobi problem. First, let us compute 
the coordinate expression of the Legendre maps associated to the Lagrangian density given by 
(l30l) . In the natural coordinates of the restricted Legendre map HC: J^vr 

J^TT* has the following coordinate expression 


Pa — 


from where the local expression of the extended Legendre map is straightforwardly deduced as 


1 


; HC p = V{x\u°‘) - - , 


Using the coordinate expression of both Legendre maps we obtain a local Hamiltonian func¬ 
tion associated to the Lagrangian density £, whose coordinate expression is 


H{x\ u'-,pf) = \ gffpip>^ - Vix\u ^), 


where, as before, 



is the inverse matrix of 



Using the Hamiltonian section h G r(/r) specified by this local Hamiltonian function we 
dehne the Hamilton-Cartan forms &h = h*Q G Qh = —dQh € J7™+^(J^7r*), whose 

coordinate expressions are 


0/j = p*„dn“ Ad™ ^Xi-{^ stfPaP^g -V{x\u°‘)] d™x , 


2 


1 




dV 


Qh = -dpi, A du" Ad™ ^Xi + ^ 1 + atfPadp^p ^ d™x 


2 du^ 


du^ 


(35) 

(36) 


Then, a representative of the class of vr^-transverse multivector fields {Xh) U j£™(J^7r*) 
which are solutions to the Hamiltonian field equation (1161) is given in coordinates by 



d 

dxl 


I ~<y.d i 

+ 9ijPa 


d 

du^ 


+ Gi 





(37) 
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with the functions j satisfying the Hamilton-De Bonder-Weyl equations (j 20 p , which in this 
case are 

m 


i=l 

in addition to the integrability conditions (f 2 T|) (if necessary). 

Let us state the generalized Hamiltonian Hamilton-Jacobi problem for this field theory. Let 
s G r( 7 rg) be a section given locally by = (x*,rt“)), where G C°°{E) are 

local functions. By Proposition [ 6 l the section s is a solution to the generalized Hamiltonian 
Hamilton-Jacobi problem if, and only if, the multivector field lYh given in coordinates by (j37p 
is tangent to the submanifold Im(s) ^ J^n* defined locally by the constraints Pa~Sa = 0 - 
Computing, the tangency condition gives rise to equations ([22]), that is, 





idg^ 





(38) 


^ + g- 
dxJ 


Pc 


M. 

' duP 




Im(s) 


= 0 . 


Now, in order to state the Hamiltonian Hamilton-Jacobi problem, we must require in addition 
the section s G r(vrg) to satisfy the condition = 0 or, equivalently, d{h o s) = 0. Bearing 
in mind the coordinate expression (I35p of the Hamilton-Cartan (m -|- l)-form the pull-back 
s*flh G is given locally by 


= 


2 du^ 


J _ 
C/3 


dV ^ ,„/3 , 


‘du^ 


dqi 

du^ A d”*x - Tv4dM^ A du" A d^-^Xi 
duP 


Hence, the condition s*B/i = 0 gives equations (f23ll for this problem, that is. 


a/3 


M + 1^ 

(9x* 2 du^ 


*3 


„i „3 
*a*/3 


dV 

du^ 


^ * 


^ = 0 

^du^ 


duP 


c) 

du'P 


Finally, we compute the Hamilton-Jacobi equation in the Hamiltonian formalism. Since 
the section s G r(7rg) is a solution to the Hamiltonian Hamilton-Jacobi problem, we have that 
hos G H^{E) is a closed form. Hence, there exists a local vr-semibasic (m—l)-form io G Q"^~^{E) 
such that du = ho s. In coordinates, the local expression of the m-form ho s is 

hos = 4dn“ A d^-^Xi - Q 44 - , 

Then, if the (m — l)-form oj is given locally by w = the condition duj = h o s gives 

rise to the following equations 


E 


dW^ 

dx^ 


ah j J 
ij *a*/3 


+ n SIT SkS' 


H(x\rt“) = 0 


dW^ _ 
dup‘ ~ 


which may be combined to give the classic Hamilton-Jacobi equation in the Hamiltonian for¬ 
malism p24p . and it coincides with the equation (|34l) obtained in the Lagrangian formulation, 
that is, 

apdW^ dW^ 

du'P duh 


^ dW^ 1 
2 = 1 


V{x\u^) = 0. 
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5.3 Minimal surfaces in dimension three 


Let us consider the following problem: we look for smooth maps </>: —>■ M such that the 

set graph((/>) C which is a surface in has minimal area and satisfies certain boundary 
conditions. This problem can be modeled as a first-order classical field theory with base manifold 
M = and 1-dimensional fibers, that is, Li = x M. Let (x, y) be the global coordinates in 
M = and {x,y,u) the global coordinates in Li = In these coordinates, the canonical 
volume form y G is given by y = dxAdy. Then, in the natural coordinates (x, y, u, ui, U 2 ) 

of J^TT, the Lagrangian density L G f?^(J^ 7 r) for this field theory is 


£ = dx A dy . 


(39) 


Observe that £ is a regular Lagrangian density, since the Hessian matrix of the Lagrangian 
function L = associated with £ and y is 

d^L \ 1 ^^2 -b 1 —UiU2^ 

" ^(ITufTuiF V-^1^2 uj + i) ’ 


duiduj 


which has determinant 


det 




duiduj J (1 -b uf -b ^ 2 )^ 


/O. 


Lagrangian formalism 

The local expression of the Cartan 2-form 0£ G l 7 ^(J^ 7 r) is 


0 ^ = 


ui 


-b uf -b ■ 


idtt A dy — 


U2 


y/l -b uf -b 


idu A dx -b 


1 


Un 


^/l + Vq + 


zdx A dy , (40) 


Un 


from where we deduce the following coordinate expression for the Cartan 3-form G J?^(J^ 7 r) 

ttlU2 


^c = - 


Uo -b 1 


\/(lTnfTu|)3 

U 1 U 2 


dtti A du A dy -b 


-b 


V^(TTufTu|)3 

Ml 


;=dMi A du A dx -b 


^(1 -bM^ -bM^)3 

nf -b 1 


dM 2 A dn A dy 


^(1 -bM^ -bM^)3 


;=dMi A dx A dy -b 


^y{l + ui + ui)^ 
m 


:^du 2 A du A dx (41) 
dM 2 A dx A dy . 


Thus, a multivector field Xc G solution to the Lagrangian field equation ([7]) is given in 

coordinates by 


Ti; = /X 1 AX 2 = / 


. f d d d ^ d , 

= f ( w b Ml— b -^1,1 w - b -^1,2w — 1 A 

Vdx du dui du 2 


f d d ^ d d 

1 w-b M 2 7 ^-b £ 2,1 w-b £ 2,2 w— 

\dy du dui du 2 


where / G C°°( J^vr) is a non-vanishing function, and the functions Fij satisfy the Euler-Lagrange 
equations for multivector fields (jlip . which in this case reduce to the following single equation 


(1 -b M 2 )£i ,1 — MiM 2 (£i ,2 + £ 2 , 1 ) + (1 + '*^i)£ 2,2 — 0 . 


(42) 


In addition, since we need the multivector fields solution to the field equation to be integrable 
in order to give a suitable Hamilton-Jacobi formulation, we must require equations ([3]) to hold. 
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In our case, since = Xi A X 2 , we must require [Xi,X 2 ] = 0. Computing, we have 

[Xi,X2 ] = (Fi,2 -F2 ,i)^ 

ou 

(dF2 1 8X2 1 ^ 8X2 1 ^ 8F2 1 

y UX UU OUi UU 2 

8 Fi 1 8 F 11 ^ 8 Fi 1 ^ 8 Fi I 

-7/,o--^2,1^--^2,2^- 


8 y ciu 


5ui 


+ 


8F2^2 8F2^2 


8 x 

8Fy2 


+ Ul' 


- U 2 


8 u 
9Fi,2 


8 F 2 2 <9-F2 2 

au 2 




ciui 

5Fi,2 


' 2 , 1 ^— 

8ui 


-F, 


2,2 


ciM 2 / 5u2 


8 y ^ 8 u 

Then, requiring this last expression to vanish, we have that the multivector fields Xc which are 
solutions to the field equation are integrable if, and only if, the following 3 equations hold 


-^ 1,2 — -^ 2,1 = 0 , 


9 ^ 2,1 , 9 ^ 2,1 , ^ 9 ^ 2,1 , ^ 5 ^ 2,1 8Fiy 

+ + Fi^i—^ + Fi^2 


8 x 

8F2,2 


+ Ul 


8 u 

8F2,2 


8ui 


8 F 22 

+ +^ 1,2 


8u2 

8F22 


8y 

8F,2 


dFi,i p 

U2— - F, 


8 x 8 u ’ 8 ui ’ 8 u 2 8 y 

A particular solution to equations (H3l) is given by 


— U2- 


8 u 

8 ^ 

8 u 


8Fy, 8F^y 

2 , 1 -?^- ^2,2^ - — U , 

OU^ OU 2 


1 


(43) 


8F\ 2 8F\ 2 

-i"2,lW^-i"2,2W7^ = 0. 

8 u^ 


F - 


F - F - - ^“2 

Ul,2 — U2,l — -F— — -TT- 

ay ox 


F - 

X2,2 — -F— ■ 

8y 


8u2 


Moreover, one can easily check that these functions are also a solution to the Euler-Lagrange 
equation (j42p . Therefore, a particular holonomic multivector field solution to the field equation 
is given in coordinates by 


„ f 8 8 8 ui 8 8 u 2 8 

\OX ou ox OUl ox 0 U 2 


8 8 8ui 8 8 u2 8 

^ 8 y^ 8 u^ 8 y 8 ui 8 y 8 u 2 


(44) 


In the following we use the particular solution given by (|44l) with / = 1 to state the Lagrangian 
Hamilton-Jacobi problem. 

First we state the generalized version of the Hamilton-Jacobi problem. Let T € r( 7 r^) be 
a jet field locally given by 'L(x,y,u) = {x,y,u,ijji{x,y,u),il! 2 ix,y,u)), with V’i )'02 £ C°°(M^) 
being local functions. Then, from Proposition [3] we know that T is a solution to the generalized 
Lagrangian Hamilton-Jacobi problem if, and only if, the Euler-Lagrange multivector field X /2 
given locally by (14411 is tangent to the submanifold Im('k) J^vr, which is locally defined by 
the constraints Uj — 'ijjj = 0, j = 1, 2. Then, from the conditions lj{Xi){uj — 'ijjj) = 0 we obtain 
the following systems of 4 partial differential equations 


/ n 

= 0 

8u 


, 811)2 

= 0 
8u 


, ^'01 

02 -^ = 0 
8u 


, i9'02 

02 ^ = 0 , 
8u 


which admits the following set of solutions 


{0i(x,?/,n) = fi{x,y),il) 2 {x,y,u) = f 2 {x,y)}, (45) 

where /i ,/2 € C°°(M^) are functions depending only on the coordinates of the base manifold 
M = R^, that is, they are constant along the fibers of the bundle R^ —>■ R^. Observe that both 
functions may be constant, and even vanish everywhere. 
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In order to obtain the equations of the Lagrangian Hamilton-Jacobi problem, we require 
in addition the jet field 'h € r( 7 r^) to satisfy the condition = 0. Bearing in mind the 

coordinate expression of the Cartan 3-form fljc € given in (j41|] . the pull-back € 

by a jet field u) = (x,y,u,'ipi,'ip 2 ) gives 


= 


\/(i + V'l + V’i)^ 


(^^2 + 1 )^ - V ' lV ’2 


d^2 d'lpi 
dx dy 


+ + 1 ) 


94^2 

dy 


du A dx A dy . 


Therefore, = 0 if, and only if, the following partial differential equation holds 


(V’2 + “ ^lV’2 


dj^ dtpi 

dx dy 


+ + 1 ) 


d^ 

dy 


= 0 . 


It is easy to check that, from the functions in the set (|45D of solutions to the generalized 
Lagrangian Hamilton-Jacobi problem, the following functions are solutions to the Lagrangian 
Hamilton-Jacobi problem: 

{■01 = /(y),V '2 = f{x)} 

with /, / G C°°{E) being functions depending only on the coordinate x and y, respectively (with 
possibly one or both of them vanishing). In addition, when = y(x,y), we do not obtain a 
closed formula for ipj in terms of the function y, but there may be functions satisfying the arising 
partial differential equation. 


Finally, we state the Hamilton-Jacobi equation in the Lagrangian formalism. Since the 
jet held T G r( 7 r^) is a solution to the Lagrangian Hamilton-Jacobi problem, we have that 
d(\J/*0£) = 0. Thus, there exists a 1-form uj G n^{E) given locally by w = W^dy — W'^dx such 
that dcu = T*0£. The pull-back of the Cartan 2-form 0£ by T gives in coordinates 

'k*0£ = — - N^idu A dy — V' 2 du A dx -t- dx A dy) . 

vi+W+W 

Hence, requiring 'k*0£ = dw, we obtain 


dW^ dW"^ _ 1 dW^ _ il)i dW"^ 

dx ^ dy -k V’l + ' du ’ du 

which may be combined in the following single equation 

dw^ dw‘^ _ I /aiT^y 

dx ^ dy y \ 9 m y \ 9 m y 


__h__ 

y/l + ’ 


(46) 


This is the Hamilton-Jacobi equation for this held theory. 


Hamiltonian formalism 


Now we state the Hamiltonian formulation of the Hamilton-Jacobi problem for this held theory. 
In the natural coordinates (x, y, M,y^,p^) of the restricted Legendre map ^ J^vr* 

associated to the Lagrangian density C given by (1391) has the following coordinate expression 


ECp^ = 


Ml 


FC*p^ = 


-|- M^ -|- 


^^2 


Mn 


\ATMf+” 


From this last expression we deduce the coordinate expression of the extended Legendre map, 
which is 


FC = 


Ml 




FCp^ = 


U2 


u‘i + M^ 


F£ p = - 


1 


•^/TTm^T^ 
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as well as the coordinate expression of the (local) inverse map J-£ ^: J^tt* —)• J^vr 




P 


; (H£-^ru2 = 


x/T^W^W 


P 


x/T^W^W 


Observe that the inverse Legendre map is not defined on the points of J^tt* satisfying + 

(p2)2 = 1. 

The local Hamiltonian function associated to the Lagrangian formulation is then given by 

H(x,y,u,p\p^) = -x/l - (pi)2 - (p 2)2 . 


Using the Hamiltonian section h € r(p) specified by this local Hamiltonian function we define 
the Hamilton-Cartan forms Oh = h*Q G = h*Q G TT*), whose coordinate 

expressions are 

Oh = p^du Ady — p^du A dx + x/l — (p^)^ — (p^)^ dx A dp , (47) 

Qh = —dp^ A du A dp + dp^ A du A dx 

H- ^ (p^dp^ A dx A dp + p^dp^ A dx A dp) . (48) 

- (pl)2 - (p2)2 


Then, a locally decomposable 2-vector field T/j G jC^(J^7r*) solution to the field equation 
i(Alh)^h = 0 is locally given by 




A 


P 


d 


yl ^ I ^2 


d 


^9x — (pl)2 — (p2p du dp^ clp2 

f d 


+ 


p" 


dy ' Y^l - (pi)2 - (p 2)2 gu~^ ^dp^ ^ dp^ j 


with the functions Gl satisfying the Hamilton-De Bonder-Weyl equations (1201) . which in this 
case reduce to the following single equation 


G\ + Gl = d. 


(49) 


Following the same procedure given in the Lagrangian formalism, an integrability condition 
must be required to this multivector field. Prom [23] we know that a particular choice of a 
locally decomposable and integrable multivector field solution to the field equation is given in 
coordinates by 


Xh = f 



p" 


d dp^ d 

+ 


_ (pl )2 _ (p 2)2 du ' dx dp^ 
’2 d dp^ d 


+ 


dp^ d 
dx dp'^ 
dp'^ d \ 


a/1 - (pi)2 - (p2)2 du dy dp^ dy dp 


P 


(50) 


+ 


+ 


As in the Lagrangian formalism, in the following we use the particular solution given by (|50p 
with / = 1 to state the Hamiltonian Hamilton-Jacobi problem. 


In order to state the generalized Hamiltonian Hamilton-Jacobi problem, let s G r(7r^3) be a 
section given in coordinates by s(x, p, u) = (x, p, u, s^(x, p, u), s^(x, p, u)), where G C°°(M^) 
are local functions. By Proposition [6l the section s is a solution to the generalized Hamiltonian 
Hamilton-Jacobi problem if, and only if, the multivector field lYh given in coordinates by (|50p 
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is tangent to the submanifold Im(s) ^ J^tt* defined locally by the constraints jj — = 0, 

j = 1, 2. Then the tangency of Xh along the submanifold Im(s) gives the following system of 4 
partial differential equations. 

ds^ ^ ds‘^ ^ 

~ ’ ’ 

ds^ ^ ds^ ^ 

This system of partial differential equations admits the following set of local solutions 

{s^(x, y, u) = f^{x, y), s^{x, y, u) = f{x, y)} , (51) 


where /^,/^ G are functions depending only on the coordinates {x,y) of the base 

manifold M = and satisfying (/^)^ + (/^)^ 1- 

Now, to obtain the equation of the Hamiltonian Hamilton-Jacobi problem, we require in 
addition that the section s € r(7r^3) satisfies the condition s*Qh = 0 or, equivalently, we require 
the 2-form hos € Q‘^{E) to be closed. In coordinates, bearing in mind the coordinate expression 
(1481) of the Hamilton-Cartan 3-form G we have 

du A dx A dy , 

from where the condition s*rLh = 0 is locally equivalent to the equation: 

ds^ ds"^ _ ^ 
dx dy 

This additional equation restricts the set of solutions (l5T|) to the following: 

{s^ = f{y),s^ = f{x)}. 


s*nh = 


ds^ ds"^ 
dx dy 


As in the Lagrangian Hamilton-Jacobi problem, /, / G C°°{E) are functions depending only on 
the coordinate x and y, respectively, and they may vanish. In addition, when s* = g{x,y), we 
do not obtain a closed formula for in terms of the function g, but there may be functions 
satisfying the arising partial differential equation. Recall that both /, / must satisfy /^ + /^ 1. 

Finally, we state the Hamilton-Jacobi equation in the Hamiltonian formalism. Since the 
section s G r(7rg3) is a solution to the Hamiltonian Hamilton-Jacobi problem, the form ho s a 
Q‘^{E) is closed. Thus, there exists a 1-form uj G Q^{E) given locally by w = IT^dy — W‘^dx 
such that do; = ho s. In coordinates, the form ho s = s*&h is given by 

ho s = s^du A dy — s^du A dx J- a/ 1 — (s^)2 — (s^)^ dx A dy . 


Thus, requiring ho s = dw, we obtain 


dW^ dW"^ , - 

+ ^ = - (Sl)2 _ (^2)2 . 


dW^ 


dx dy 

which may be combined to obtain equation 


du 
that is. 


= s 


dW^ 

du 


= s 


dW^ dW"^ 


dx 


+ 


dy 


fm 
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6 Conclusions and further research 

Starting from the geometric Hamilton-Jacobi theory developed mainly in mm for mechanical 
systems and using the results given in [39] as standpoint, we have stated a geometric framework 
for first-order classical field theories described in the multisymplectic setting. 

The theory has been developed for the Lagrangian and the Hamiltonian formalisms. In 
both cases, first we have stated the so-called generalized Hamilton-Jacobi problem, which is the 
most natural one in this geometrical ambient, and hence we have defined from it the standard 
Hamilton-Jacobi problem. Particular solutions to these problems are dehned and characterized 
in several equivalent ways and, in particular, one of these characterizations for the standard case 
in the Hamiltonian formalism, when written in natural coordinates, leads the classical Hamilton- 
Jacobi equation for field theories. After that, the definition and geometric characterization of 
complete solutions is also given and different features about them are discussed. Finally, the 
equivalence between the Lagrangian and the Hamiltonian Hamilton-Jacobi problems is also 
proved. 

It is important to point out that this generalization of the theory has been achieved using 
distributions in the jet bundles and multimomentum bundles where the Lagrangian and Hamil¬ 
tonian formalisms of multisymplectic classical field theories are developed. These are integrable 
distributions whose integral sections are the solutions to the Lagrangian and Hamiltonian field 
equations, and they are represented, in general, by means of equivalence classes of multivector 
fields. This choice has enabled us to give a construction of the Hamilton-Jacobi theory in a very 
natural way. Thus, our model is different from that given by L. Vitagliano in |59] for higher- 
order field theories, who uses connections as the main geometrical tool and a unified formalism 
to describe the Lagrangian and Hamiltonian formulations at once. 

We have analyzed several examples. First, non-autonomous mechanical systems can be 
considered as a special situation of held theories, and hence the Hamilton-Jacobi equation for 
these systems has been recovered from our model as a particular case. Second, we have applied 
our results to obtain the Hamilton-Jacobi equation for held theories described by quadratic affine 
Lagrangians. Finally, we have written this equation for a more particular example: minimal 
surfaces in dimension three. 

As further research, we believe that our geometric framework for the Hamilton-Jacobi theory 
can be extended to higher-order held theories using the formulations in [9] [52| , thus generalizing 
the results of [181119] for higher-order mechanics and giving a different but equivalent perspective 
to that of [59] for this kind of theories. 

A very relevant application of the Hamilton-Jacobi theory for hrst-order held theories would 
be to the Palatini approach of General Relativity and, once the extension to higher-order held 
theory is made, also to the Einstein-Hilbert Lagrangian approach, as well as to other gravita¬ 
tional theories. As a previous step, a suitable multisymplectic description of these gravitational 
models must be done and, although there are some recent attempts to do it [58] . more work in 
this way is necessary and research in this way is in progress. 

Another interesting question in the ambient of the Hamilton-Jacobi theory is the existence 
of conserved quantities (or conservation laws), and the integrability of the system. It is known 
(see [13] for more details) that, in the case of mechanics, for a dynamical system with n degrees 
of freedom, the existence of complete solutions to the generalized Hamilton-Jacobi problem is 
associated with the local existence of families of n functions which are constants of motion. From 
a geometrical perspective, a complete solution is a foliation in the fiber bundle which represents 
the phase space of the system, which is transverse to the fibers, and such that the dynamical 
vector field is tangent to the leaves of this foliation. Then, these leaves are locally the level 
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sets of the functions which are constants of motion. In addition, a complete solution to the 
Hamilton-Jacobi problem corresponds to a Lagrangian foliation (with respect to the symplectic 
structure, canonical or not, which the phase space is endowed with), and thus the constants of 
motion are in involution and the system is completely integrable. 

In the case of field theories, from a geometrical point of view, the situation is quite similar: 
as it is defined in Sections 13.31 and 14.31 complete solutions to the generalized Hamilton-Jacobi 
problem endow the jet and multimomentum bundles with foliations which are transverse to the 
fibers and such that their leaves contain the image of the sections solution to the field equations. 
When we consider just the Hamilton-Jacobi problem, these are m-Lagrangian foliations (with 
respect to the corresponding multisymplectic structures, in the sense defined in [12]). Never¬ 
theless, up to our knowledge, the notion of “integrability” is not clearly stated in these cases. 
Furthermore, although the leaves of these foliations can be also locally defined as level sets of 
families of functions, how to associate these functions with conservation laws in field theories 
must be investigated and, even if the foliation is m-Lagrangian, these functions are not said 
“to be in involution” because, although there are several attempts to define unambiguously a 
Poisson bracket (for functions) in covariant field theories (that is, in multisymplectic geometry), 
this problem is not solved in a completely satisfactory way. The discussion on all these topics is 
also under research. 
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